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Abstract 

We study the electrical susceptibility of a hydrogen gas at equilibrium, partially ionized 
by thermal excitations. The gas is described as a quantum plasma of point protons and 
electrons, interacting via the Coulomb potential. Using the newly developped diagram- 
matical technique of screened cluster expansions, we calculate exactly the wavenumber- 
dependent susceptibility in the atomic limit, where most charges are bound into hydrogen 
atoms. A transition from conductive to dielectric behaviour occurs when the wave length 
is decreased well below the Debye screening length. The standard formula for the di- 
electric function of an ideal gas of hydrogen atoms is recovered in an appropriate scaling 
limit. The derivation treats all effects arising from the Coulomb interaction (screening, 
binding, polarization) in a fully coherent way, without intermediate approximation nor 
modelization. 

PACS: 52.25.Mq; 52.25.Jm; 

Keywords: Quantum plasma, atomic limit, dielectric function, screened cluster expan- 
sion 



1 Introduction 

Screening in a non-relativistic nucleo-electronic plasma is one of the most important con- 
sequence of the long range of the Coulomb force. For quantum mechanical point charges, 
there is a number of physical effects that can participate to the screening mechanisms. 
First, there is a collective screening phenomenon which leads to the formation of neutral- 
izing polarization clouds (as in the classical plasma). The clouds are made of unbound 
(ionized) charges, and are always present at any non zero temperature. The spatial ex- 
tension of such clouds is of the order, say, of the Debye-Hiickel screening length. On the 
scale of the Bohr radius, chemical binding may lead to the formation of neutral atoms 
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or molecules. Finally, in a partially recombined plasma, there exists also a dielectric 
screening due to the polarization of atomic dipoles and molecules. 

In usual theories, these different phenomena lead to different ways of modelling the 
system, each of them having its own range of validity. In a conducting phase, the ionic 
screening due to the unbound charges may be taken into account by replacing in calcu- 
lations the bare Coulomb potential by a mean field potential obtained within the Debye- 
Hiickel or Random Phase Approximation. If one is interested in a dielectric phase, it is 
often appropriate to neglect ionic screening and to separate the problem into two parts. 
The role of quantum mechanics is limited to an a priori calculation of the polarizability 
of a single atom, and then the problem is treated in the framework of classical statistical 
mechanics of preformed dipoles, characterized by these quantum mechanical atomic data. 
It is nevertheless true that all these effects have a single common origin, the Coulomb in- 
teraction. All of them should stem in a consistent way from the basic A^-body Hamiltonian 

^ Ip.p ^ 1 

describing the Coulomb interaction of N quantum particles (point nuclei and electrons) 
of species a = 1,2, .. . with charges Cq, and masses rria- Such a fundamental attitude is 
legitimate not only as a matter of principle, but also because of the need to a have a 
coherent scheme for understanding the interplay and relative importance of these effects. 

In this paper, we present a detailed study of the response function x(r) of a quantum 
plasma to a classical localized external charge density Cext(r). Its Fourier transform 

^ few 

is defined as the ratio, to linear order, of the induced charge Cind(k) in the plasma to the 
external charge Ccxt(^) at wave length k. It is related to the dielectric function e(k) by 

X(fc)=e-i(fe)-l (3) 

so that purely metallic behaviour (e(0) = oo) is characterized by the perfect screening 
relation 

lim x(fc) = -1 (4) 

The response function can be expressed in terms of the charge fluctuations in the plasma 
by 



X(fc) = -^ iyrS{k,T). (5) 



W\ 



In 0, S{k,T) is the imaginary time displaced charge charge correlation function (the 
charge structure factor of the quantum plasma at imaginary time r), and /? = (/c^T)"^ 
is the inverse temperature. These formulae hold in a fluid phase and it will be assumed 
throughout this paper that the system is spatially uniform. The relation @ together 
with is the quantum analogue of the Stillinger-Lovett perfect screening condition for 
a classical plasma p. 

In a mean field treatment, the low wave number behaviour of x(^) can be represented 
by the Debye-Hiickel type formula 
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where k~ = Ad is the Debye length. In an uniform state of the quantum system, it 
turns out that the relation @ is always true at any non-zero temperature. This result 
can be obtained from an analysis of the constraints imposed by the hierarchy equations 
for imaginary time Green's functions [2 1^ • It can as well be derived in the formalism of 
charged loops (section 4) using various schemes for summing Mayer graphs j^. Hence, 
an infinitely extended quantum plasma is formally always conducting, even in a phase 
composed mainly of neutral entities (atoms, molecules). This is due to the tiny amount 
of free charges that are present by thermal ionization. 

The main question we adress in this paper is: under what conditions does the system 
of quantum charges exhibit a dielectric behaviour ? To keep the discussion reasonnably 
simple, we restrict it to the electron-proton (e-p) system with A'p protons, A'e electrons 
and Hamiltonian H]\f^]\[^. If electrons and protons recombine into a dilute gas of hydrogen 
atoms with density p^, one expects, at an elementary level, the system to be characterized 
by a dielectric constant 

e ~ 1 + A-KpnOH (7) 

where Oh is the polarizability of the hydrogen atom. In order to establish such a result 
starting from the many-body Hamiltonian Hn^j^^, it is necessary to give a precise meaning 
to the recombination of protons and electrons into hydrogen atoms. This is formulated in 
the so-called atomic limit described in section 2. In this limit one lets the temperature and 
the density tend to zero in a coupled way. Low temperature favors binding over ionization, 
whereas low density, by increasing the available phase space, favors dissociation. The rate 
at which the density is reduced as T — > determines an entropy-energy balance that 
selects the formation of certain chemical species. If this rate is within a certain range, 
hydrogen atoms are formed, and it can rigorously be shown that the equation of state 
becomes asymptotic to that of an ideal gas of hydrogen atoms El • 

In order to put the issues in proper perspective, we give in section 3 a naive but 
mathematically ill defined derivation of ((7|, disregarding all collective screening effects 
tied with the long range of the Coulomb potential. This also enables us to specify the 
range of wave numbers k for which ^ is expected to hold: in view of (jlj) k should not be 
too small to avoid the perfect screening regime, but also not too large so that the atom 
experiences an uniform perturbation on the scale of the extension of its center of mass 
wave packet (see (|l9|) V 

The rest of the paper is devoted to an exact derivation of Q in the atomic limit. The 
wavenumber-dependent response function x(fc) does not appear to be easily analyzed with 
the rigorous methods used in [SI [7]. Here, we use the technique of quantum Mayer graphs 
that gives a straightforward expansion for two particle correlation functions needed to 
calculate x(fc). It allows a derivation of Q that is exact in the sense that it does not 
involve any intermediate model or approximation. In particular the many-body problem 
is fully taken into account and the existence of preformed atoms is not assumed. However 
the derivation remains formal to the extent that results are established for each individ- 
ual graph, without control of the convergence of the diagrammatic series. The present 
technique is also suited to explore the vicinity of the atomic limit and to systematically 
calculate corrections to the ideal gas behaviour. For instance, non-ideal contributions to 
the Saha equation of state are derived in 

In section 4 we briefly recall the loop representation of the Coulomb gas. This formal- 
ism arises from the Feynman-Kac path integral representation of the quantum mechanical 
Gibbs factor. Charges of the same type and same statistics are collected according to a 
permutation cycle of length q {q = 1,2, . . .) into a random path (a Brownian path) called 
a g-loop (for a review and references, see In terms of loops, the statistical averages 
are performed according to the same rules as in classical statistical mechanics. As a con- 
sequence the powerful method of Mayer graphs is available in the space of loops. Thus 
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an ensemble of point quantum mechanical charges becomes isomorphic to a classical-like 
system of fluctuating multipoles. Two points can then conveniently be made at this stage. 
First, the divergences due to the long range of the Coulomb potential can be cured by 
the introduction, via partial resummations, of an effective screened potential (which is 
the quantum analogue of the usual Debye potential). The properties of this potential are 
studied in details in jlUj . Morover the formalism offers the possibility of an easy derivation 
of the response function since the rules of classical linear response apply. In particular, 
in the langage of charged loops, the perfect screening relation (jlj) takes the same form as 
the Stillinger-Lovett rule for a classical system of structured ions (see H81|) ) . 

The loop expansion needs to be properly reorganized to perform estimations in the 
atomic limit. The reason is that an element of the space of loops consists into a number 
of charges of the same species, therefore not capable to bind. By a diagrammatical 
reorganization, the loop expansion can be converted into the screened cluster expansion 
(section 5). The details of this reorganization are presented in In the latter expansion, 
the basic elements are all possible clusters of positive and negative charges, thus candidates 
for atomic or molecular recombination. Let us just mention that for a system of particles 
interacting via a short-range potential, this expansion exactly coincides with the usual 
quantum mechanical virial expansion. In the Coulomb case, the expansion undergoes 
a number of modifications because of the necessity to deal with an effective screened 
potential. 

At this point, we have in hand the necessary tools to investigate the response function 
in the atomic limit. In section 6, we select and study the leading graphs yielding the 
dielectric response. There are essentially two aspects to be controlled. As the density goes 
to zero, the screened potential reduces to the (non integrable) bare Coulomb potential. 
One must make sure that this does not create any divergence in the atomic limit. Secondly, 
one must control that all effects due to excited states and ionized states of the hydrogen 
atom become negligible when the temperature vanishes. We find it worth to give complete 
proofs of these non trivial mathematical points (technical parts of the proofs are relegated 
in appendix A). We check then that all the other graphs involving clusters of more than 
two particles do not contribute in the limit. For this we rely on the general analysis 
presented in . Eventually, section 7 is devoted to a discussion of the response function 
in the wave number region that interpolates between perfect and dielectric screening. 
Our results are compared with the extended RPA dielectric function introduced in the 
framework of standard many-body perturbation theory. 



2 The atomic limit 

The notion of atom or molecule in the many-body problem can only receive a precise 
meaning in an asymptotic sense. The temperature must tend to zero to give a predominant 
weight to bound states over ionized states and the density p should be small enough 
(p-i/3 ^ ^]^g Bohr radius), to have spatially non overlapping atoms. To formulate the 
atomic limit, we introduce first the grand canonical densities of three independent species, 
the electrons (e), the protons (p) and the hydrogen atoms (a) in their ground state. 



Pf = ( o "^fe2 ) exp(-/3(^j -//j)), j = e, p, a (8) 



3/2 

\2Trph? ) 

where \ij are the respective chemical potentials, ra.^ = nie + rup, = Ep = 0, Ei^ < 
is the ground state energy of the hydrogen atom, and dj is the spin degeneracy (da = 4, 
de = dp = 2). Here all the effects of the Coulomb interaction are disregarded except for 
binding energy \Eg^\ of the atom. 
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The law of chemical equilibrium for the dissociation reaction e + p a requires 



Hi, = tJ-e + fJ-p (9) 
and one also must have charge neutrality = p^. Introducing the combinations 

/^ = ^— , ^ = ^— (10) 

it is easily seen that the neutrality condition imposes the choice 

u = m = j^i^^ = o{p-'). (11) 

Hence (jHJ becomes 
with 



K = hJ^ (14) 

V "^a 

the thermal de Broglie length of a particle of mass m^. The ideal density of an atom, 
molecule or ion with A^p protons and A'e electrons and ground state energy Ej^^p^^ is of 
the form 

Pn^n. = ,^ t2"''\3/2 exp[-/3(i?iVp7V. - ^^{Ne + iVp) + 0{p-'))] (15) 

where we used the fact that fieNg + fipNp = fi{Ne + A'^p) + ^'{Ne — Np) and (|llj) has been 
taken into account. XNpNc is the corresponding thermal length and dj^pNc degeneracy 
factor. 

The situation where the atomic density dominates all other ionic or molecular densities 
at low temperature is characterized by 



pL^»pj^p^., f^^oc. (16) 



More generally, we require that the probability of occurrence of hydrogen atoms dominates 
that of all possible other configurations of protons and electrons. This will happen if there 
exists a range I of chemical potential ^ such that all the inequalities 

< - 2/x < SjVpJVe - KNe + Np), (N,, Np) / (0, 0), (1, 1) (17) 

can be simultaneously satisfied, where E^^^^ is the infimum of the spectrum of H^^n^- 
The range / (if any) is located above E^. since setting {Nc,Np) = (1,0) (single electron) 
in (fT7|) gives fi > £'a. Moreover setting {Nc,Np) = (2,2) (hydrogen molecule) implies 
< \{E22 — Eq). Note that E'a < \{E22 — Eg) since the binding energy [£'22 — 2£'a| gained 
by the formation of a hydrogen molecule is less than the binding energy of the hydrogen 
atom himself, a well known fact. These two cases do not exhaust all the constraints 
imposed by the inequalities (|17|) . From now on we make the plausible, but yet unproven, 
hypothesis that there exists an interval /a =]£'a5 -E'a + A], A > 0, such that ((T7|) holds 
when /X S /a- 
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It is worthwhile to note that the vaUdity of this hypothesis is equivalent to the possi- 
bility to find an optimal constant for the stability of matter estimate that we write here 
in the form 

HN^N. >-B{N, + Np-l), (iVe,A^p)/(0,0),(l,l) (18) 

for some positive constant B. It is conjectured (but not proven) that there exists a 
stability constant B strictly less than |£'a| for all cases except of course for the hydrogen 
atom itself. It is then easily checked that the hypothesis B < \Ea.\, {N^, Np) ^ (0, 0), (1, 1) 
in H18|) is equivalent to the possibility of finding the non void interval /a such that all the 
inequalities (|17)1 are satisfied, /i G /a; hence assuring the existence of the atomic phase. 
Under this condition, it has been proven that El E] 

piP,^i) = pl^{l + 0{e-'^)), /iG/A, /3^oo (19) 

namely that the grand canonical pressure p{P, n) verifies the equation of state of a perfect 
gas of hydrogen atoms of density p'^ up to an exponentially small correction in /?. This 
is the situation we will consider in this paper and is referred to as the atomic limit 
/i G /a, P ^ oo. 

In general, collective screening effects are provided not only by free electrons and 
protons, but also by all types of more complex ions that have appreciable densities PjVpAfe 
for given values of the temperature and the chemical potential. In this paper, we study 
the interplay between dielectric and ionic screening in the situation when the latter is 
predominantly due to unbound electrons and protons. This will certainly be the case if p 
is chosen sufficiently close to £'a, in a possibly smaller interval Is =]£'a,-E'a + S] C Ia so 
that 

0<E^-2p<-p< En^n, - p{Np + iVe) pe Is (20) 

holds for all {Np,Nc) ^ (0, 0), (1, 0), (0, 1), (1, 1). We consider hence an atomic phase 
under the conditions that the most probable entities in the system are, after hydrogen 
atoms, ionized electrons and protons, namely 

pL' » Pf = Pp » diN., (Np, iVe) ^ (0, 0), (1, 0), (0, 1), (1, 1). (21) 

In the quantum Mayer graphs expansion, the collective screening effects are embodied into 
the effective screened potential obtained by chain resummations (see section 5). Under 
the conditions just described, we take as intermediate chain points those corresponding 
to individual electrons and protons. Then the resulting screening length in © will 
only involve the densities of free charges H12|). If it turns out, in an other regime, that 
collective screening effects mainly arise from other types of ions, another definition of the 
effective potential must be adopted in consequence. 



3 Elementary description of screening 

In this section, we describe, in an elementary (and non rigorous) manner, the screening of a 
classical external charge by an e-p plasma in the atomic limit. This heuristic description 
is intended to serve as a gentle introduction to the calculation of the next sections, in 
which the response function is studied in full generality. 

In the atomic limit, the e-p plasma approaches an ideal gas of hydrogen atoms of den- 
sity Pa^. Under the influence of an external charge, the hydrogen atoms get polarized and 
screen partially the external charge. At low temperatures and low densities, one expects 
this effect to be described, at leading order, by the dielectric constant 0. According 
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to © and at linear order in p^^, the system's response function should therefore approach 



X{k) ~ -^^p'ia^ (22) 

in an appropriate range of wave numbers. 

In view of Q, let us make a straightforward low activity expansion of the imaginary 
time displaced charge charge correlation function S{k,s). This correlation is defined by 

S{k,s) = {c{k,s)c{0)) (23) 

where (. . .) is the grand canonical average. For a particle system evolving with the 
Hamiltonian Hat 

N 

CN{r) = Y,eaAr-r^) (24) 

i=l 

is the microscopic charge density operator, and 

N 

CN{k,s) = e^H^( J^e„,e-**''->)e-^«^ (25) 

i=l 

is its Fourier transform at "imaginary time" s. In general, if A^v are A^-particle observ- 
ables, the first terms of the low fugacity expansion of the grand-canonical average of A, 
with 



(A) = ^ E ' - = E ^Tre-Z^H- (26) 

N=l ' N=0 



are 



(A) = zTrAie-^^i + ^ixr Age'^^^ _ (Tr^^e-'^^i)(Tr e"^"!)] + 0{z^) (27) 
In the e-p system, the one-particle and two-particles Hamiltonians are 

p2 

H„ = a = e,p (28) 

2 2 

Haiaa = 7j^ + 7j^ + e«ie«2^(ri -ra) (29) 



It is understood in (|27)) that the traces include a summation on the particle species and 
are carried out on properly antisymmetrized electronic and protonic quantum states. 

We apply ^ to ^ and are interested in extracting from it the contribution pro- 
portional to p!f (see which comprises the factor exp[— /3£'a] where -Ea is the ground 
state energy of the hydrogen atom. Such a contribution must come from the second term 
of (|27|) . the only one involving the two particle Hamiltonian. This contribution is then 

\ E E Tre-^""i"2C2(fc,s)c2(0) (30) 

where we have made the sum over species and spin explicit and the trace is now purely 
configurational (zq, = exp[/3/ia] is the fugacity of the particles of species a). Moreover, 
only the part of H3U|) pertaining to an electron-proton pair with Hamiltonian Hep will 
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involve the factor exp[— /3£'a]. So keeping only the terms ai / Q2 in (|Hn|l gives with (j211), 
H25() and /ic + /^p = 2/Lt 

4e2/3/^e2Tr{e-('3-")^^p(e-*'"-'"p - e-^*'-^^)e-^^^p(5(rp) - (5(re))}. (31) 

We now proceed to the evaluation of this trace, going to the center of mass variables. 
Since 



Hep = :^ + :^-eV(re-rp) (32) 

separates into a center of mass Hamiltonian Hcm = P^/2M and a relative Hamiltonian 
H = p^/2m — e^/|r|, its eigenstates are of the form \p,r\) = \p) ^ |n), with eigenvalue 
Ep + En , where 

Hcm Ip) = Ep \p) and H | n) = | n) (33) 

(n is a generic index for the eigenstates (including the ionized states), n = being the 
ground state). This allows one to write the e-p contribution to S{k,s) as 

Ae^Pf^e^ Up I dp' e-(^-'^)[^''+^"le-^[^p'+^n'l(p,n|e-*'=-'-i'-e-''=-'^|p',n'> 

•' n,n'>0 

X {p\n'\5{vp) - 5{r,)\p,n) (34) 



where rp = R — and = ~^ ^'^^ the quantum operators associated to the 
position of the proton and of the electron. The eigenstates of Hcm are of course plane 
waves {R\p) = (27r)^^/^ exp[ip • R\ with energy Ep = h^p'^/2M. The matrix elements 
in (|34|) are therefore equal to 



(p, n| e-^'=-^p - e-*'^-^^ \p' , n'> = 5{p' -p- k)A„n'{k) (35) 

1 

(2V)= 

where we have defined 



{p' = p + fe, n'l 5(rp) - 5(re) |p, n) = -——A„.„{-k) (36) 



A„w{k) = (n| e^*'-^'" - e'^^M^ |n'> . (37) 
From the above matrix elements, the e-p contribution to x{^) becomes 

=2 rP 



Xep(fc) = -^4e2/^'^ / dr/,/^(fc) e-('3-^)^"e--^n'|An'(fc)l 

n,n'>0 



(38) 



where 



(39) 



Uk) ^J^l dpe-^(i-^)^-e-^^^^+^ 
The Gaussian integral in fs{k) can be calculated with the result 

where Aa is the thermal wave length of the atom. In the atomic limit, /3 is very large 
and the dominant terms in (|38|) are those with n = or n' = 0, because they contain 
the exponentially growing factor exp[—(3EQ], which is much greater than exp[— /?i?n] if 
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n > (£^0 = -E'a)- Keeping only these ground state terms and factoring out the atomic 
density we find that the hydrogen atoms give the dominant contribution 

Xep(fc) ^ -^T^f^ia^{k, P), f3^oo, (41) 

to the response function, where 

a,(k,l3) = ^ Td^e-^'^'^'^^^-^) V e-'^(^"-^«)e-^(^n'-^")|An'(fc)p. (42) 

nn'=0 

The quantity aH(fc,/3) can be interpreted as the polarizabihty of a hydrogen atom at 
inverse temperature /3 ^ {Ei —Eq)~^ in an external electric field varying on the scale k~^. 
Let us discuss briefly the result (|42j) for aii{k, (3), before commenting on the method used 
to derive it. 

When A: — > 0, the external perturbation becomes uniform on larger and larger dis- 
tances, so that one expects to recover the polarizabihty of a hydrogen atom in an 
uniform electric field. For this, it is necessary to have the perturbation uniform on the 
scale of the dispersion of the center of mass distribution, namely k <^ A~^. If it is the case, 
the factor exp[— ^A;^A^^(1 — ^)] in 1)421) can be approximated by one and the r-integral 
in (|42j) is easily evaluated: 

2e2 ^ 1 - e-f3{En-Eo) 
a^[kJ)^-^Y. ^ ^ \Mn{k)\\ fcAa«l. (43) 

In obtaining (|4,'-{|) . we used the symmetry properties of the matrix element: A„„r{k) = 
^n'n(~^) ^-ncl |ylnn/(fc)p = |Ann'(— fc)P (this follows from the fact that H commutes with 
the parity operation). Moreover, for /? large, Aa is much greater than the Bohr radius, 
and hence k ^ A^^^ ^ a^^, so that one can make the dipole approximation 

\Ao„{k)\^ ^k^\{0\k-r\n)\'^, k <^ a-\ (44) 

Notice that (0| r |0) =0 by parity. For k <C A^;^ and /3 > (E^i - Eq)'^, auik^/S) reduces 
therefore indeed to the ground state polarizabihty of a hydrogen atom in an uniform 
electric field [T^ : 



(0| fc • r |n 



2 
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aH(fc,/3) ^ an = = -a|, k (45) 



n>l 



(The term n = in (|l3|) is of order j3e^{ka-s)^ /k"^ = a%{kf5e^){ka^) and can be neglected 
when k <^ A^ ^.) This will give the anticipated form (|22() of the response function provided 
that k is not in a range where perfect screening is prevailing. In view of ©, we must 
require 



K 



2 



« 4np^au, (46) 



A;2 + K 
or, equivalently. 



fe^»^|.A-. (47) 



In gni), 



K = ./87r/3eVe (48) 
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is the usual inverse Debye length for an e-p plasma and Eq. (|47|) defines a length Ai which 
is the border-line scale below which ionic screening starts taking place. 

Summarizing the discussion, we expect that in the vicinity of the atomic limit, the re- 
sponse function of the e-p system should be essentially constant and equal to the value 
in the range 

XY^'^k<.X-\ (49) 

This region corresponds to the dielectric behaviour of the e-p system in its atomic phase. 
As k is further decreased, ionic screening will predominate, and interpolates between 
the values —Anp^^au and —1. Notice finally that ^ as /c ^ oo since atoms cannot 
be polarized under the effect of an electric field with infinitely fast spatial oscillations. 

The cross-over length Ai between dielectric and ionic screening was derived in an 
infinitely extended state of the e-p system. The same length can also be obtained by 
considering the screening of a point external charge cq immersed in a spherical sample of 
radius L. We assume L much smaller than the Debye screening length (L <^ k^^). The 
charge density of free charges induced around eo is given in the Debye-Hiickel approxima- 
tion by Cfi.ee(^) = — cqK'^ exp(— Kr)/(47rr). The screening provided by the free charges is 
hence 

[ dr Cfreeir) = -eo [ dune"" ~ -eo^(KL)^. (50) 

J\r\<L Jo ^ 

Comparing with the induced charge due to polarisation — 47rpiJ^aHeO) the sample is indeed 
expected to display a dielectric behaviour if <^ Sirp'^au/K'^ = Aj, in agreement with 
the cross-over distance defined in H47() . 

The present elementary description of screening shows that the dielectric regime of 
the response function is formally characterized by 

x(^) 

lim — rr > —Airan, kX^ —>■ 0, A;Ai oo, fx ^ I. (51) 

This should be contrasted with a strict zero temperature limit at fixed k. If one takes the 
limit /3 ^ oo of an{k,/3) at fixed k (let s = t/P in (|^ and note that ^fc^A^ = PE^), a 
Dirac delta function appears: 

lim (]e-'^iE.-^o+s{E„,-E.-E^)-E^s^] ^ 25(^„ -Eo + s{E„, - - E^) - Eks"^) 

/3—>cx} 

= — Ero('5) + — ^^i^-^)' n • n = 0. (52) 

Ek + En' — Eq Ek + En — Eq 

Using ds5{s) = 1/2 and the symmetry properties of the matrix element |^on(^)|5 we 
find 

„2 I /i,m2 



2e^^ \Aon{k)\ 

n>0 

If we now let A; ^ in (|53() , we obtain 



lim QH(fc,/3) = °" ^ = aH(fc,oo). (53) 



hm aH(fc, oo) = aH + ^^^^3 ' ^^1 |0) ? (54) 

which is obviously different from (|45|) . The result (|54|) can also be obtained from a direct 
calculation in the ground state of Hgp (|32|) . The additional contribution in H54j) has its 
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origin in the coupling of the center of mass of the atom with the external electric field. 
This contribution does not vanish in the limit /c — > 0, because the external field does not 
appear to be uniform (even when k is very small) to the atoms, the latter being entirely 
delocalized at zero temperature (here k~^ <C Aa = oo). 

We comment now on the above calculation of xik)'- in fact all steps of this calculation 
are ill-defined, except for the final result (|41|) and (|42|) . First of all, since x(^) represents 
the response to an external charge density localized in the bulk, the infinite volume limit 
should be taken first to disregard all boundary effects. But for a non-integrable potential, 
all virial coefficients of the series (|27|) (except the first one) diverge in the thermodynamic 
limit. The first term in (|27() gives the contribution (the factor 2 accounts for the spin 
degeneracy) 

-^2 Ve^z^ f'ds (0| e^^'^e-^-e-^'^ |0) = (55) 
to the response function, where we defined the function 

Ak)^4.p^el-^^ rd.e-^^^^^^(i-). (56) 

„ i^T^KY' Jo 

Notice that K^{k) reduces at A: = to (see (fl8|) l. using the definitions ^ of the ideal 
densities and the neutrality condition pjj^ = p]^. Hence, when kXa ^ 1, K'^{k) ~ k^, and 
(|55j) is essentially identical to © as long as k ^ k. However, since (|55|) diverges when k 
0, while x(^) goes to —1 in this limit, there must obviously exist some other contributions, 
coming from n-body states in the series (|27() . which cannot be neglected even at low 
density. Physically, these many-body contributions are associated to a collective effect in 
the system: the screening of the Coulomb interaction by "screening clouds" of particles 
appearing around the charges. This effect must be taken into account systematically in 
order to cure the problem of the Coulomb divergences in ((77|) . 

Furthermore, besides the collective screening effects, we have not dealt with the more 
complex entities (e.g. the hydrogen molecule) that can be formed by chemical binding. In 
view of (|16|) and (|17j) . such entities should not contribute to the dielectric response in the 
regime defined by the condition ()49() in the atomic limit. Finally, one should check that 
in this limit, excited and ionized states of the hydrogen atom do not contribute either. 
The main goal of the following sections is to provide a calculation algorithm which is free 
from these difficulties. 



4 The loop representation of x(^) 

We recall in this section the loop representation of the response function x(^)- That 
representation is obtained when the grand partition function Ha of the quantum plasma 
is written in a classical form, by using the Feynman-Kac path integral formula J3] and 
collecting permutations with the same cycle structures [TH [T^ : 

N=0 ' 1=1 

This so-called magic formula relies on the following definitions. The element of phase 
space £, called a loop, is a collection of q particles of the same species exchanged in a 
cycle. A loop 

£ = {R, a, q,X{s)), 0<s<q (58) 
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is specified by its position R in space, a particle species a, a number of particles q, and a 
shape X{s) with X{0) = X{q) = 0. The closed path 



R{s) = R + XaX{s) (59) 

describes the trajectories of the q particles, which are located at R{k — l),k = l,...,q 
as they move from their position to the position of the next particle in the cycle in unit 
"time". The path X(s) is distributed according to a normalized Gaussian measure T>{X), 
with covariance 

D{X)X^{s)X,{t) = 6^,,q\mm(-,-) - (60) 

L \q qJ q'^i 

Integration over phase space means integration over space and summation over all internal 
degrees of freedom of the loop (which we denote collectively by x = ['^iQiX): 

d£---= I ArI dx---= / ARY.fl I D(^)--- (61) 

J J J a=l q=l 



The interaction energy of loops is the sum of two-body interaction potentials U (£i, Cn) 

i<j 



"^i^i^j with the interaction between two different loops 



V{Ci,Cj) = ea,eo,^ r ds dt 6is - t)V{Ri{s) - Rjit)). (62) 
Jo Jo 

In (inH), 6{s) = S{s — n) is the Dirac comb of period one. V{Ci, Cj) is hence the 

sum of the Coulomb interactions between the particles in the loop jOi and the particles in 
the loop Cj as they move along their trajectory. The loop potential is clearly a function 
of the relative distance Ri — Rj and of the internal constitution of the loops: 

V{CuCj) = V{Ri-Rj,Xi,X,)- (63) 
Eventually, the activity of a loop reads 



(the factor 2 takes the spin degeneracy into account) where 

U{J^) = 4 r dsi r ds2 (1 - - S2)V{R{si) - R{S2)) (65) 

^ Jo Jo 

is the sum of the mutual interactions of the particles within a loop (the factor (1 — 
(5[sj] [sj]) excludes the self-energies of the q particles). The above rules define the statistical 
mechanics of the system of charged loops, which we call the loop representation of the 
quantum plasma. Note that the interaction potential (|62|) inherited from the Feynman- 
Kac formula is not equal to the electrostatic interaction between two classical charged 
wires, which would be 

VeieciCi, Lj) = e«,e«^. ^ ds dt V{Ri{s) - Rj{t)) (66) 
Jo Jo 

Although the formalism of loops has a classical structure, the difference between V{Ci, Cj) 
(|62j) and VeiedCi, Cj) is responsible for the absence of exponential screening in the quan- 
tum plasma jSj. This difference is the occurrence of the equal time condition 5{s — t) 
which characterizes the quantum mechanical aspect of the interaction (|62() . 
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In the loop representation, one can define the loop correlation functions according to 
the usual definitions. Introducing the loop density p{C) = the average loop 

density and the two-loop distribution function are 

p{C) = {p{C)) , P.{Ci,C2) = (/5(£i)p(£2)) - m2)) , (67) 

where the average is taken with respect to the statistical ensemble of loops defined in (|57|) . 
and coincident points are included in pT^[Ca, Cf,). It is appropriate to recall that the charge 
sum rule holds in the system of loops [n 



y dr y dxaQaeaa PT{r,Xa,Xb) =0. (68) 

Any fixed loop of charge Qbea,, is surrounded by a cloud of loops of opposite total average 
charge. 

When a localized external charge Cext(*') is immersed in the plasma, the corresponding 
induced charge density Ci^dif) can be calculated, within the formalism of loops, according 
to the rules of classical linear response theory. The external potential 

Kxt(r-) = / dr'^^^ (69) 



due to Cext('') is represented in the system of loops by 

VeM= rdsVeAR + XaX{s)). (70) 
Jo 

It gives rise to an interaction energy 

C/ext = J d£e„p(£)Fext(>C) (71) 

which has to be added to the pair interaction of loops U in (|57() . Then the response of 
the loop density pind(^) in presence of the external potential l|70j) is given at linear order 
by the standard formula 

Pind(-C) = -/3 y ea'V,,t{C')p^{C, C) (72) 

To form the induced charge density at r, one has to integrate p\nd{C-) on the shape of the 
loop and take into account that a loop carries total charge Caq: 

Cind(r-) = -/3 j dR' J dx J dx'e^qe^, J^' dsVe.tiR' + Xa'X'{s))p^{r,x,R',x')- (73) 

The final expression for the dimensionless response function ((2} follows after Fourier trans- 
formation of (|7,''{|) ■ taking (|()9|) into account, 

X{^) = -^ j ^XaJ dxbea^ea, f^' dne'^-^^''^^^^^Mk,Xa,Xb). (74) 

Here pTik, Xa, Xb) is the Fourier transform of the translation invariant truncated loop-loop 
density fluctuation PT{f,Xa,Xb)- The formula H74() has flrst been derived by Cornu P^ . 

There is an interesting interpretation of the perfect screening relation @ in terms of 
the statistical mechanics of random charged loops. We split xi^) into two parts 

H^) = -^(^'°°P(fe) + M'°°P(fe)) (75) 
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where 



5''°°P(fc) = y dxa j dxbeaaqaeatQbPT{k,Xa,Xb) (76) 

is the Fourier transform of the charge-charge correlation of loops and 

Mi°op(fc) = JdxaJ dxbea^Qaea, J^'" dr(e^'=-^'^'>^''(-) - 1) M^, Xa, Xb) (77) 

comprises the multipolar contributions of the loops to x(^)- Because of the charge 
sum rule (|68)1 together with rotational invariance, both S^°°P(fc) and M'°°P(fc) have to 
be 0(|fcp) as k — > 0. The rotational symmetry forces these terms to take the form 

5'°°P(fc) ~ Ifcl^i j dr |r|25'°°P(r), A; ^ (78) 



and 



M>°°P(fc) ~ i J dXadXbea^Qaik ■ d{xb)){k ■ VkPT{k,Xa,Xb)\k=o) 



drr-pi°°P(r) = / dr|r|2V-P^°°P(r), A: ^ 0. (79) 



6 

In ((7^. we have defined the dipole of a loop by d{x) = eaq)^a Jq dr X{t) and introduced 
the polarization vector 

pi°°P(r) = [dxa [ dxbea^qad{xb)pAr,Xa,Xb)- (80) 



as equal to the charge-dipole correlation of loops. With these definitions, the perfect 
screening relation written in terms of charge-charge and charge-dipole correlation of loops 
takes the classical form of the second moment Stillinger-Lovett condition, namely 



~6~ 



J dr |r|2(5i°°P(r) - V • P'°°P(r)) = -1. (81) 



The same relation holds in classical models of structured ions where both the charge 
density and the polarization charge —V • P{r) participate to the constitution of the 
screening cloud 



5 The screened cluster expansion 

The loop formalism leads itself naturally to the introduction of Mayer graphs on the space 
of loops. A vertex receives the weight z{C) (|BH) and a bond the factor exp[—pV{Ci,Cj)] — l. 
Since the loop pair potential (|62j) behaves as the Coulomb potential qiea^qjeaJ\Ri — Rj\ 
itself, the bonds are not integrable at large distances, and partial resummations are needed. 
All divergencies can be removed by introducing the effective screened potential (t)[Ca-,Cb) = 
4>{Ra — Rb, Xa, Xb) defined as the sum of chains 

-P(p{Ci,Cj)= O O + O • O + O • • O + •■• (82) 

where the bond is the linearized Mayer bond —f3V{Ci, Cj). This potential is the quantum 
analogue of the classical Debye potential (see [TOi;). For the purpose of this paper, it is 
convenient to use a slightly different definition of (j) by restricting the intermediate (black) 
points to one-particle loops. Its properties are the same as those of the effective potential 
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studied in JO]- The potential cp describes the screening effects due to ionized protons and 
electrons with the inverse Debye screening length k H46() . At short distances \r\ <^ k~^, (j) 
reduces to the bare Coulomb potential V{Ca,^b) between loops. At distances |r| ~ k~^, 
(j) ~ gaCQ^gbeajj exp[— Kr]/r approaches the standard Debye potential that describes the 
classical collective screening effects. At large distances, (j) has an tail corresponding 
to dipole-dipole interactions between the loops (this tail is responsible for the algebraic 
decays of the correlations in the quantum plasma). 

Once all chain summations (jSH) are performed in the Mayer graphs, the resulting 
prototype graphs still obey Mayer diagrammatic rules, with two kinds of bonds {Fc = —/?</> 
and Fr, = exp [—/?(/)] — 1 + two kinds of weights, and two additional rules that prevent 
double counting (see HO]). Like in a system of classical dipoles, the bond Fc{Ca,J~-b) is 
at the borderline of integrability. The prototype graphs are integrable at large distances 
provided that the integrations on the internal variables of the loops (the shape X{s)) are 
performed first. 

Though the prototype graphs describe non perturbatively screened Coulomb interac- 
tions between quantum particles, they are not adapted to the evaluation of the response 
function x(^) i^i the atomic limit. Indeed, the small parameters in this limit are the ideal 
densities (section ISJ 

pid^^^ oc e-^(^^p^e-^(^p+^=)) « 1 (83) 

which are not apparent (or even present) in individual prototype graphs. This problem 
occurs because the vertices in prototype graphs involve single loops, which are groups 
of particles of the same species, with only repulsive pair-wise interactions (included in 
the loop-activity ()64() '). Attractive interactions appear in the graphs as bonds connecting 
electronic and protonic loops, but a given graph does not in general include the total set 
of pair-wise interactions that would be necessary to form bound states between electrons 
and protons (see JJ). Individual loop-Mayer graphs are therefore not in direct correspon- 
dence with the ideal atomic or molecular densities occurring in the atomic limit. For this 
reason, it is convenient to collect prototype graphs to form new graphs involving clusters 
of protons and electrons, together with all their mutual interactions and proper statistics, 
in such a way that the new effective activities approach the ideal densities in the atomic 
limit. This reorganization, called the screened cluster expansion, is worked out in details 
in We recall here only the final diagrammatic rules, taking into account the minor 
modifications introduced by the choice of a slightly different effective potential (p. We 
stress that this expansion is nothing but the usual quantum cluster expansion, suitably 
generalized to take into account the screening due to the long range of the Coulomb po- 
tential. 

The screened cluster expansion for equilibrium quantities of the quantum e-p plasma 
is expressed in terms of Mayer graphs G, with the following definitions of bonds and 
weights ^n\ . 

Vertices A vertex C in a graph G is a cluster of particles, denoted by C{Np,Nc,), 
containing A'^p protons and A^e electrons. The internal state of a cluster involves all 
possible partitions of the protons and electrons into sets of protonic and electronic loops. 
Let 

La 

Qa = [qi,---,qLj, ^qi = N^ (84) 

i=l 

be a partition of Na into La subsets of qk particles, k = 1, . . . ,La, with qi > q2 > ■ ■ ■ > 
qi^. Here La runs from 1 to Na- To a partition ((5p,(5e) of the Ap protons and Ae 
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electrons, we associate a cluster of loops 

C(Qp, Qe) = {4'^ , 4'^ 41} (85) 

where carries q"^^ particles of species a (k = 1, L^). The variables associated to a 
cluster C{Np,Nc) of a graph G are Qpj Qe, ^(Qp, Qe)- The statistical weight of a cluster 
reads 

where na{q) is the number of loops containing g particles of species a in the partition Qa- 
In (|HH|l. 2:</,(-C) is a renormalized loop activity 

z4£) = z(£)e^«(^), /r(£) = ^(F - 0)(£,£), (87) 

and the truncated Mayer coefficient B'^ ^ is defined by a suitable truncation of the usual 
Mayer coefficient B(f,^N for loops with pair interactions (p (see jllj. This truncation 
ensures that B'^ ^ remains integrable over the relative distances between the loops when 
(j) is replaced by V . The two first truncated Mayer coefficients are B^ i = 1 and 

Bl^ = exp(-/3<A) + + (88) 

A vertex corresponding to a cluster C where the associated variables are not integrated 
over is called a root point (or white point). The integration over an internal (or black) 
point is performed according to the measure 

^(^) = E / fl ^^t^ ft I n D(^?^) ft D(^?) (89) 

Q^n^'' k=l k=l '' k=l k=l 



Bonds Two clusters Cj and Cj are connected by at most one bond J-(j,{Ci,Cj) which 
can be either /3^<I>^/2! or -/3^<I>^/3!. The potential $(Ci, Cj) is the total interaction 

potential between the loop clusters Ci{Q^f^ jQ^f') and Cj{Q^^\Q^^^) that describe the 
internal states of Cj and Cj respectively, i.e. 

HC^,C,) = m,Cj} = E E '^(^'^')- (90) 



Special rules In two cases, the weight (|86j) of a black cluster must be modified to 
avoid double counting: 

(i) If C is an intermediate cluster in a convolution (— /3^>) * (— /3<^) and contains only a 
single electron or proton (represented by a loop C with q = 1), then its weight is 

Z;{C) = z^C) - z{C) (91) 

instead of Z'^{C) = z<^(£). 

(ii) If C is a cluster connected to the rest of the graph by a sing le bond l{(3^f and 
contains only a single electron or proton, then its weight is also given by (|91|1. 
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Diagrammatic expansion of the two-body loop density In view of (|7I|) . we 

need the screened cluster expansion of pT^{Ca,Ch). According to eq. (4.31) of it is 
given by 

X I l[B{Ck)z;iCk)[l[:F^]^ 

k 

+ D(Ca)D(a)( Y ^i^^^^a) ^^^k.C,)) 

X Zl{Ca)Zl{C,) j\{j^{Ct)Zl{Ck)[\[T^]^ (92) 

k 

where the graphs G involve either a single root cluster Cab that contains both loops 
Ca and Ch (which can be the same loop since coincident points are included), or two 
roots clusters Ca and with Ca in Ca and in C^. The symmetry factor Sg is the 
number of permutations of labelled black clusters that leave the product of bonds [f]^ ^<j>]G 
unchanged (only clusters with identical numbers of protons and electrons are permuted). 
The sum runs over all topologically different unlabelled graphs C which are no longer 
integrable over the relative distances between the clusters {Cj}i=i^,,,^n when (/> is replaced 
by V . Except for this additional constraint (represented by the star in X^g), the graphs 
G have the same topological structure as the familiar Mayer diagrams: the one-particle 
points are replaced by particle clusters and the usual Mayer links are now the bonds J^^. 
A few graphs occuring in the screened cluster expansion of x(^) are shown on Fig. 1. 




Figure 1: Examples of graphs occuring in the screened cluster expansion of x(fc). The clusters are drawn as 
circles containing a certain number of protons and electrons. The graphs have either one root (white) cluster 

Cab (graphs a, b, c and e), or two root clusters Ca and Cb (with Ca drawn to the left of Cf,). The bonds , 

= and ^ are -(3^, (/3<i>)^/2 and (-/3$)^/3! respectively. 



Because of the central role played by the effective potential 1)82(1 . it is useful to recall 
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from the following exact formula for its Fourier transform 
^(fc, Xa, Xb) = y exp[-ik ■ r] 4>{r, Xa, Xb) 

= e^^e^, Tel., rd.,e^'=-[^«^»(-)-^^^^(^^)l V 

Jo Jo + f^^ik,n 



^-i2TTn{sa-si,) 



(93) 

where the screening factor for frequency n is 



(94) 



This is the expression (33) of jTHj restricted to one-particle loops. The functional integral 
can be evaluated using H6U() with the result 

^Hk,n) = A.pY.^l-^^^ /'d.e-^^^^^^(i-)e^2™. (95) 

^ {^T^Kr Jo 

Notice that the zero frequency term n^ik^n = 0) = ^^{k) is identical to 1)56^ . 

6 Dielectric screening in the atomic limit 

The screened cluster expansion allows one to study the response function of the e-p plasma 
in the atomic limit, without encountering any divergence, thanks to the screening effects 
embodied in (f). We show here, by a term by term analysis of the diagrammatic series, that 
the response function tends in the coupled limit (3 oo, k ^ {fj, £ Is) defined in to 
the value —ATrp^^Un, as expected for a gas of hydrogen atoms at low density. This result, 
announced in constitutes a first principles derivation of a dielectric constant in the 
quantum e-p plasma. 

In the atomic limit, the fugacities Za = exp[/?/XQ,] (and hence the densities p'^) vanish 
exponentially fast as /3 ^ cxd. The different lengths in the system are therefore ordered 
according to 

A„ <C /3e2 < (p!f )-i/3 < (96) 

where the first inequality follows from oc y/p. In the present limit, the screening length 
diverges and the effective potential (j) tends to the bare Coulomb potential V. The 
difference 4>ch = (j) — V between cj) and V is therefore expected to be small. It is indeed 
proven in that at sufficiently small fugacities 



\<Pc\i{r,Xa,Xb)\ < KqaQbe^K, (97) 

where X is a constant independent of the loop variables. This allows us to replace in all 
graphs the renormalized loop activities Z(j,{C) (|H7jl by the bare activities z{C), since 

e^H(^) = l + 0{f3e'^K). (98) 

Similarly, it will be legitimate at leading order to replace (j)hj V in all statistical weights 

z;{c). 
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6.1 The mean field ionic contribution 



We start by calculating the simplest contributions to xik), namely that of free protons 
and electrons (graphs (a) and (b) of Fig. 1). The statistical weight 1)86(1 of the root cluster 
Cab = C(1,0) or C(0, 1) involves a single protonic or electronic loop C = {R,a,q = 1,^) 
with activity z^{C) z{C) = 2za/ (27rA^)^/2. According to (fTl]) and (jHU, the contribution 
to x(fc) of these graphs is 



M V / Df£) C dse'^-^-^'^'^ 



Eq. (|99() is identical to the first term (|55() of the "nai've" virial expansion. The mean field 
(or Debye-Hiickel) result Q is obtained by adding to (j99|) the contributions of the four 
graphs 







+ 



+ 







(100) 



At leading order, (|100|) becomes 



V e„„e,, / D(£JD(£,) /' dTe^^-^''«''(^)z(£„)z(£5)(-/3</'(fc, Xa, X^))- (101) 

„=e.r> J Jo 



aa=e,p 
ot,=e,p 



Using (P^ . and 

Tdr /"^dse*2™ /" D(£)z(£)e*'=-^[^M-^(^)] =5^,0 /"^ dr /" B{^) z{£)e'''-^^^^\ (102) 
Jo Jo J ' Jo J 

which follows from the periodicity of £(s) and of exp[i27rnr], the contribution (|101(1 is 
easily evaluated, and gives, when added to (|99|) . 

XMF(fc) = + fc2[fc2 + ^2(fc)] - -fc2 + ^2(fc)- (lO^) 

This result agrees with the Debye-Hiickel formula © since «;2(/i;) — > in the limit 
(|51j) . Notice that ()103j) saturates the perfect screening relation Q . The sum of all other 
contributions to x(^) must therefore vanish at A: = 0. 



6.2 The atomic contribution 

The dielectric screening effect ()41|) is expected to be contained in the graph consisting 
in a single root cluster C(l, 1) (graph (c) of Fig. 1), since the latter graph describes an 
interacting electron-proton pair (H atom). We calculate here the contribution of this 
graph, with a precise mathematical control on its asymptotic value in the atomic limit. 
The cluster C(l,l) is made of one protonic loop Cp = (rp,p, l,£p) and one electronic 
loop Ce = (T*e,e, l,£e)- ^^s Contribution reads 



47r/3 



dr^l drp I D(£J I D(£p) e,e, j\s e^^'^M^h-'^-^^ 5{r ,) 

X z^iCe)z^{Cp)l3l^{Ce,Cp). (104) 



i=e,p 
i=e,P 
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Using the translation invariance B'^ 2(^0, ^p) = 2(^6 — fp,Xe,Xp) a-^d introducing the 
change of variables Tp — > — r, we can rewrite ()1U4|) as 



In the atomic limit, the factors exp[— go to 1, and the effective potential (f), which 
enters the Mayer coefficient i3J 21 tends to the bare Coulomb potential V{r,Xa,Xb)- Recall 
that the truncation in B"^ 2 such that i3J 2 remains integrable when (p is replaced by V. 
The integral in (|105|) is hence finite for any k, but since it is divided by k'^ and A; ^ 
in the dielectric limit defined in it is important to note that the integral behaves at 
small k as k'^ times 

Jdr j D(0 I D(^p)(fc.er).(fc.e[r + AeC(5)-Ap^p(s)])^J,2(r,Xc,Xp). (106) 

Since the integrand in the above expression becomes non integrable at zero density (it 
decays as l/r"^ when cj) is replaced by V), we further truncate B'^2 by introducing the 
decomposition 

Bl2{^a, A) = B^liCa, Cb) + ^ i-m^a, A))' + | {-P'PiCa, Ck)f , (107) 

which defines ;SJ'^(£a, £b). Now B'^\{Ca-,Cb) decays as \Ra — -Rfcl"^, so that second 
moments are finite. The terms cf)'^ and cff' give a negligible contribution in the atomic 
limit, as compared to the free charges contribution (see section 6.3). After a little algebra, 
the difference between B'^2 a-^^d ^Sy^ is found to be 

BlliCa, A) = Bll{Ca, A) + R1 + R2 + R3, (108) 

where the remaining terms Ri are 

R,{Ca,Cb) = (e-'^'^ch(^.,A) _ i)Bl7^^^Ca,Cb) (109) 

R2iCa,Cb)= V ^ — i-PMCa,Cf,)n-PViCa,Cb)r (HO) 



nl m\ 

n,m=0 
n+m>6 



CO ^ 5 
R3{Ca,Cb) = y^-A-PMCa,Ct)r y^{-PV{Ca,Ch)r. (Ill 

n=6 m=0 

We define Xa.t{k) as (|1()5|) with bare activities and with B'^2 replaced by By^. 



gifc-Ao4<,(s) _j_ g«fc-Ap€p(s) _ g-ifc-rgifc Ap4p(s) _ gifc-rgifc-Ac€c(s) \ _ (112) 



X 

To evaluate (|1()5|) in the atomic limit, we must establish the following two points: 

(i) Replacing cj) by V. We must show that (|112j) is indeed the leading behaviour of p()5|) 
in the atomic limit. This is done in appendix (section A.l), where we prove that the 
contribution of the remaining terms Ri in (|1()8|) are o(/9^), and hence negligible as 
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(ii) Excited and ionized states. We must show that the contributions to x(fc) of ex- 
cited and ionized states of the hydrogen atom, which are included in H112() . are also 
negligible in the atomic limit. This is done in section A. 2 of the appendix. 

Using point (i), we proceed to the evaluation of Xat(^) at leading order. We introduce 
the relative and center of mass coordinates, defined by 

(x,^^{s) = X,Us) + ^Xm ^^3^ 

lAp^p(s) = Aa^a(s)-^A^(s) 

with 

"^=^T^' M = me + mp (114) 
mm 

Notice that AeAp = AaA. It is easy to verify that the Gaussian measure D(^g)D(^p) and 
D(^j^)D(^) have the same covariance. In these new variables, the functional integrations 
over and ^ factorise. The center of mass integration 
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(2,A2)3/2 / my"--'^'^^'^ = fsik) (115) 
gives the factor fs{k) already encountered in (|l0|) . Xat(fc) hence becomes 

Xat(fc) = -^4e2^^^'d./,(fc) JdrJ m)Bll{r,XO 
where, according to (|1U7() . we have 

Bl^lir, A^) = exp [(3e^ C dsV{r+Xi{s))\ 1 ((3e^ f dsV{r + X^{s))\ > 0. (117) 



In order to determine to low temperature limit of (|116|) in terms of atomic eigenvalues 
and eigenstates, we convert this expression back into operator's language. Notice that the 
factor in brackets in (|116() can be rewritten as 

Defining the operator A = e* — e~' as in (|37j) . we find from the Feynman-Kac 
formula ^31 that (|ll()j) is equivalent to 

Xat(fc) = -^4e2/^'^ dr/,/^(fe)Tr{u(/3-r)AtU(r)A 

rd^^(.))"]}. (119) 



_e-/3Ho^ 

n=0 

In H119|l . the trace runs over the spectrum of the hydrogen Hamiltonian 

2 2 

H = Ho + V, Ho = |-, V = -^, V = -V, (120) 

2m |q| 

U(s) = exp[— sH] is the evolution operator, and the freely time-evolved operators are 
defined by 

A(s) = e""«Ae-""o, V(s) = e"^oVe-"^o. (121) 
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The subtraction of the freely evolving quantities in ensures the finiteness of the 

trace. The dominant low temperature terms will come from the ground state contribution 
of U(s) when evaluating the trace. Let P = |0) (0| be the projector on the ground state, 
Q = H — P = X]m>il'^) ("^1' decompose U(s) = Up(s) + Uq(s) accordingly, i.e. 
Up(s) = exp[-£;os]P and Uq(s) = U(s)Q = QU(s). We split (fnojl into 

Xat(fc)=xit^(fc) + xii^(fc), (122) 

where 

= -^4e2^^^''dT/,/^(fc)Tv{Up(/3-r)AtUp(r)A 

+ Uq(/3 - r)AtUp(r)A + Up(/3 - t)A^Vq{t)a} (123) 

is the part of Xa.t{k) which has contributions from the ground state, and xit^(^) involves 
only contributions from excited and ionized states of the hydrogen atom, as well as the 
truncations terms in (|119|) . It is shown in appendix (section A. 2) that the part xit^(^) is 
negligible in the atomic limit as compared to Xa?(fc)- In terms of atomic eigenvalues and 
eigenstates, H123() becomes 

= fdrf^/pik) Y: e-(^--)^"e-^"'|A.„,(/^)P (124) 

^ n,n'>0 
n-n'=0 

where An^'ik) is defined in Factoring out the atomic density (fT^ . we find {Eq = Ei) 



xS^(fe) = -47rpLdaH(fe,/3), (125) 

where a^{k,j3) ~ an if ^ ^ (see (jl^ V The graph Cab = C'(l)l) does therefore 
indeed describe, at leading order in the atomic limit, the dielectric screening due to the 
polarisation of the hydrogen atoms, in agreement with the anticipated result (|41|) . We 
stress that no divergences occur in the present derivation of (|125|1 , in contradistinction to 
the elementary calculation of section 121 

6.3 Higher order contributions 

We consider all graphs in the screened cluster expansion different from the pure atomic 
graph Cab = C{\,1), and argue that they give higher order contributions to x(^) in 
the atomic limit [3 oo, fi E 1^, if k 3> X^^ (see (jS^)- In fact, the estimations of 
those graphs differ from the estimations presented in 8 for the particles densities pa = 
J2g Q j D(X)p(£) only by the different weight attached to the root points. The definitions 
of weights for black clusters and bonds are indeed the same, so that the factors arising 
from integrations over these clusters can be evaluated in the same way. In |Hj , an analysis 
of the behaviour of Coulomb multiparticle Green functions shows that these factors vanish 
exponentially fast as /3 — > oo. We do not reproduce this mathematical analysis here, but 
merely adapt the arguments to the present situation. 

We assume without loss of generality that the set of fugacities {ze,Zp} satisfy the 
pseudo-neutrality condition 
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so that electrical neutrality p^^ = p^p of the ideal system (no Coulomb interactions) holds. 
With this choice, the difference between the chemical potentials pc and pp are given by 
and depends on the fugacity 

^ = e^^ p=^^^. (127) 

The choice of fugacities satisfying the pseudo-neutrality condition will allow to greatly 
reduce the number of graphs contributing at a given order. Moreover, we choose p > Eg^ 
in the interval Is so that, according to (|2H) . collective screening effects are due to ionized 
electrons and protons with densities p^^ = p^^ (see end of section I^J . The estimations 
will be performed in terms of the exponentially small factor p^^ (x z = exp[/3^], p ^ Is, 
disregarding any power law dependence in /?. The condition (PT|l together with k » Aj^^, 
or equivalently 

P ^ ^2 ^ "sPa > l^^»J 

will allow simple estimations of the graphs that avoid considering their small k behaviour. 

(a) The graph consists in the single root cluster C{Np,Ne) 

We assume (A^p, A'c) 7^ (1, 0), (0, 1), (1, 1), since these cases have already been consid- 
ered previously. The contribution of the root cluster is 

47r/? f v-^ 




A;2 



nr. r-r^ -^0 



= ^9n,nX^), (129) 

At leading order, the effective potential (/> in the renormalized activities and in the 
Mayer coefficients IS^^ can be replaced by the bare potential Y . The loop integrals in 
QNpNaik) do indeed converge at zero density (for any A; > 0), despite the long range of the 
Coulomb interaction, because of the truncation built in By. With these replacements, all 
density dependences in gNpNc (^) are contained in the prefactors of the integrals, which are 
obviously of order z^p'^^". To evaluate the low temperature behaviour of the integrals, 
we convert them back in operator's language using the Feynman-Kac formula. Similarly 
to H119() . QNpNcik) is given by a trace evaluated over suitably antisymmetrized states of the 
A'^p protons and A'^e electrons (because of the sum over the partitions), of a time ordered 
product of the Gibbs operator exp[—(3ILN^Nc] a-^d operators exp[ik ■ qj(s)], where qj(s) is 
the time evolved position operator for the i^^ particle. Notice that the truncation terms 
ensure the finiteness of this trace despite the long range of the Coulomb interaction. A 
typical term in the truncation involves Gibbs operators for sub-clusters of (Mp,Me) 7^ 
{Np,Ne) particles, Mp < A^p, Mg < N^. As /3 ^ 00, the leading behaviour of the 
truncated trace is controlled by the ground state contribution of exp[— /JHAr^jVc], which is 
proportional to exp[—f3Eiyp]y^] (discarding powers of /3). As we have shown in details in 
the case of the atomic contribution (see appendix A), the excited states and the truncation 
terms are expected to give exponentially smaller contributions at low temperatures. The 
function gNpN^if^) is hence expected to behave as 

ffJVpJVe(fc) oce-^(^^p^'=-^(^p+^^», P^oo, (130) 

where E^^j^^ is the infimum of the spectrum of KjypNc- Using (PT|) . gNpN^ik) is therefore 
bounded for (3 large by an expression exponentially smaller than the density of ionized 
charges: 

9NpnM < Pfe'^^ , (A^p,Afe)/(l,0),(0,l),(l,l), (131) 
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where F is a positive constant. In the low temperature behaviour of the truncated 
traces of Gibbs operators exp[— /^H^VpATc] are studied in details, and shown to satisfy the 
upper bound (|131|) . This bound, combined with l|128jl . is sufficient to show that the graphs 
consisting of a single root cluster C{Np,N(,) / C(l,l) do not contribute to the response 
function in the dielectric regime characterized by the limit (|51|) . 



(b) A black cluster Cio 
cluster C(A^p, N^) 



or Coi is connected by a bond — /?$ to a root 





+ 





(132) 



The integration over the black cluster involves the expression 
^ y"di?y"D(X)z(£)e^«(^)(-MAA)) 



(133) 



a=e,p 



where C = (a, q = 1,R, X) is the loop in the black cluster and £j is one of the loops in the 
cluster CNpNc- Eq. H133|) can be evaluated by using the Fourier transform (p{p,XiXi) 
with wave-numbers |p| 0. Using the low p behaviour K?{p,n) = n^Snfi + 7nP^ + C(p^) 
and the rotational invariance of z exp[I/j(£)], we find that the contributions of the terms 
n 7^ vanish by parity in the limit p — > 0: 



lim / D(X)ze^«('^)^ 



As / ds' (e*P-^"^(") 

JQ 



An 



1) 



+ 7„ 



i27rn(s— s') 



0. (134) 



In (|134|) . the two subtractions —1 could be freely introduced because n ^ 0. Only the 
term n = does therefore contribute to H133|l . and we find 



(PHI) 



■E 

Q=e,p 



(27rA2)3/2 



D(X)e^«(^) = OiPe^K) 



(135) 



The estimate 0{f3e'^K) is obtained by using and and noting that the term 

of order 1 vanishes because of the pseudo- neutrality condition (|126j) . The two graphs 
H132|) give hence, according to H135|) . a contribution exponentially smaller than the one 
associated to the graph (|129j) . Notice that without the choice of pseudo- neutrality, (|135|) 
would be of order one, and "dressing" points in a graph with a black cluster Cio or Cqi 
connected by a bond — /3$ would not raise the order in density. 

The leading behaviour (|135|) can be obtained more directly by replacing in (|133j) the 
effective loop potential (p{Ci,C2) by the Debye potential 



qiq2- 



exp[— K|r|] 



Ri — R2 



(136) 



The result H135|) follows then from the integral J dr exp[— Kr]/r = Att/k^. The replacement 
of (j) by (|136j) is a priori valid only for distances r ~ (see section but it can be 
used at all distances to evaluate the leading behaviour of (|133|) (see ^J). Indeed, at 
short distances r < A, (p{Ci,C2) is given at lowest order by the bare Coulomb potential 
V{Ci,C2) (see (jHZI))- The contribution of the region r < A to the integral (|133|1 is of 
order A^, as can be seen after the change of variables r = Xx. Since A <C k~^, this 
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contribution is small as compared to the contribution of the region A < r < k~ \ which 
is of order At very large distances, r ^ k'^, (j) decays as a dipolar potential. 

This term vanishes however by parity after integration over X, and the following term 
in the multipolar expansion is proportional to A^/r^. This term is integrable and the 
contributions to the integral (|l.S3j) of the large distances r ^ are hence also negligible 
(they are of order kA^), as the short-range contributions. It is therefore legitimate to 
replace (/> by (|136|) at all distances, since we correctly capture in this way the leading 
contributions to the integral, which come from the intermediate regions r ~ We 
stress that the exponentially growing factor 1/ k? in the estimate (|135() has its origin in 
the fact that the bond becomes non integrable as the density goes to zero in the atomic 
limit ((/)—> y). Because of this density dependence of the bonds in the graphs, the lowest 
order at which a graph contributes cannot be deduced immediately from the number of 
black clusters it contains. 



(c) A black cluster C(l, 1) is connected by a bond 

C(iVp,A^e): 





to the root cluster 



(137) 



For an electron-proton pair forming a hydrogen atom H136|l should be replaced by 

— K|r| g— K|r+a| - 



r + a 



(138) 



where a is of the order of the Bohr radius. The bond (|138)) remains integrable as k ^ 
due to neutrality and rotational invariance (the average dipole moment ea of the atom 
is zero). The graph gives hence a contribution to x(^) that is exponentially smaller than 
(|129|) because of the additional factor pif arising from the weight of the black cluster. 



(d) A black cluster C(Mp,Me) is connected by a bond — /9$ to a root 
cluster C(A^p,iVe): 




(139) 



The leading contributions associated to the integral over the center of mass of the black 
cluster arise as before from inter-cluster distances of the order of where the loop 
potential can be replaced by ()136|) . At these distances, the bond —(3^ is insensitive to the 
precise form of the loops within the clusters, since Aq <C The integration over the 

center of mass provides as before a factor jSe^ /k?, while the integrations over the relative 
distances in the cluster yield at low temperatures a factor p^VpAfe' '^^ich arise from the 
ground state contribution in the associated truncated trace. This provides the estimate 
f3e^ p^^^jq^l for the integration over the black cluster in (|139|) . Using (|2T|) . we conclude 
that the integration over the black cluster C^pN^ gives, if (Mp,Me) ^ (0, 1), (1,0), (1, 1), 
a factor vanishing exponentially fast as /3 — > oo. Notice that this estimate is also valid, 
but not optimal, in the case of a neutral cluster (see point (c) above). 

(e) The cases (b), (c) and (d) when the bond is |(/?$)^ or |(-/3$)^, for 
instance: 




(140) 
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We can argue that these graphs give smaller contributions than those discussed in (b), 
(c) and (d). This is due to the fact that these bonds give less divergent integrals as 
K ^ than the bond — /3^>. The integration of a bond provides for example a 

factor proportional to k,~^ oc z~^^'^. This can be seen by replacing as before the effective 
potential by the Debye potential (|136() . and evaluating 




(141) 



In the case of the bond we have to use the correct behaviour of (p at short 

distances (r < A) to estimate the integral of the bond, because the Debye form ()136() 
would lead to logarithmic divergences. As shown in the result is that the integral 
grows like In(KA) as /9 — > oo. The factor arising from the integration over a cluster 
connected by a bond or |(— are therefore smaller than the corresponding 

factor for a bond — and the graphs considered in (b), (c) and (d) are negligible in 
the atomic limit. With the decomposition (|107j) of iSJj; we introduced two graphs with 
bonds ||(— or ^{—Pcj))^ between an electron and a proton. Since these bonds are 
integrable at zero density^, such graphs are of order jl? (discarding powers of /?), and 
are hence negligible in the dielectric limit H51|l . in view of ()128() . 

The above analysis can be generalized to estimate the order of magnitude of an arbi- 
trary graph with one root cluster and black clusters connected by bonds — ^(/3$)^ or 
g(— (see JJ). We stress that the neutrality of the hydrogen atom plays an impor- 
tant role in the estimations (see cases (a) and (c)). In particular, when the root cluster 
C{Np,Nc) is the electron-proton cluster C(l,l), a rough estimation of its contribution 
based on H13U|) would give /Je^pjf //c^, which is not small compared to /?!^ when k is in 
the dielectric range (namely (|128j) no longer applies). But in fact, as shown in details in 
section IH?^ the root cluster C(l, 1) is really of the order /ojf Oh in the dielectric regime: the 
dangerous factor l/Zc^ is cancelled by neutrality. More generally, the estimate H13U|) is not 
optimal when C{Np, N^) is a neutral root cluster (A^'p = Ac) corresponding to a molecular 
state of Ap protons and Ae electrons. Its contribution does not behave as jie^ p^^^^^/k"^ , as 

estimated in (a), but rather as P^MpNc^NpNe with aNpN^ the polarisability of the molecule. 
If C(A^p, Ae) is not neutral, we cannot improve on the estimate f3e^ p^^^^^/k"^ and hence 
the use of (PT|) cannot be avoided when pJvpAfe density of a charged ion. 

Graphs with two root clusters We consider eventually the case of graphs con- 
taining two root clusters Ca and Ch- Since we perform in (|92|) an integration over Ca, 
the latter cluster can be treated as a black cluster with a special weight Z^*(Ca). Writ- 
ing Caa^a,, = e^sgn(eQ,^)sgn(eQ,J in (|74j) . all graphs obviously have a prefactor [3e^ /k"^, 
as before. The special weight Z^*(Ca) differs from Z'^(Ca) by an additionnal factor 
ex.p[—ik ■ Ri]sgn{eai)qi for the loop identified to Ca- These factors do not depend on f3, 
and change the estimations only in one case: the pseudo-neutrality condition cannot be 
used anymore, because of the factor sgn(ea), to cancel the contributions of order as- 
sociated to clusters Ca(l, 0) or Ca{0, 1) connected by a single bond — by adding these 
two contributions. However, because of the factor exp[—ik ■ R], the integration over the 
cluster Ca involves the Fourier transform of the bond which is of order 

/3e^pi^/{k^ + ^). (142) 

When k is large enough for (|128() to hold, the graphs with two root clusters are hence all 
negligible in the atomic limit. 

^They decrease as and at large distances and the Coulomb singularity at the origin is smoothed out 
by the functional integration — see 
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In conclusion, we have shown that if k = kp remains in the range Xj <^ k <^ \ 
the only graph giving a non vanishing contribution to x(^) / P^a limit /? ^ oo is the 

pure atomic graph C(l, 1) calculated in section 

7 On the cross-over regime 

When k is not large in front of Aj"^, ionic screening cannot be neglected, and xi^) iii" 
terpolates between the values x(0) = — 1 and —Anp^^au- This cross-over regime between 
ionic and dielectric screening can be described approximately by adding further diagrams 
to the graphs calculated in the previous section. We noted earlier that the four graphs 
which give the mean-field result saturate the perfect screening relation (jlj). 

Taken individually, the other graphs in the screened cluster expansion do not give however 
contributions to that vanish when A: — > 0. Contributions compatible with perfect 
screening are obtained by "dressing" the root points with a cluster C(l, 0) or C(0, 1) con- 
nected by a bond —P^, similarly to the "dressing" procedure used to prove the screening 
sum rules obeyed by the particle correlations using their Mayer-loop expansions |lUj . 

We can obtain a model that describes the cross-over regime between ionic and dielec- 
tric screening, by adding to xmf(^) the contributions of the "dressed" hydrogen cluster 
C(l,l): 




(143) 



where the unspecified particle in the root cluster Ca and Cb can be either a proton or an 
electron. Let us evaluate at leading order the second graph in p43p . for k <C A^^. Its 
contribution can be represented by the eight diagrams^ 




(144) 



where the vertices represent loops £ = {R,a,q = 1,^) of weight z{C) = 22;Q,/(27rA^)'^/^. 
The species of each loop (a = e or p) is indicated on the graphs. The wiggly line represents 
the bond and the root loop Cb ^ind the internal loop are connected by the bond 

every graph of (|144|1 . the integration over Ca (see (|74|) ) involves the evaluation 
of the Fourier transform of the bond —/?</>: 

D(^J I di?„e-**=-^(-M£„,£,)) (145) 

^The representation (|144|l is in fact nothing but the set of the loop-Mayer graphs in the expansion of 
PT{£a,Cb) that have been collected together in ,11| to form the considered screened cluster graph. 
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where Ci is the white loop Ch or the black loop. At leading order, (|145j) can be calculated 
by retaining only the contribution of the term n = in H93|)^. with the result: 



The contribution of the eight graphs H144() to x(^) is hence 

fc2 (27rA2)3/2 (2^A2)3/2 + ^^k) J J ^^^^^^^^P^ 

ds dt 2(^, Xe, Xp) (^e^'=-^^^-(")e-^'=-^^^-(*) + e^'=-^p^p(*)e-*'^-^p^p(*) 

_ g-«'«-^g*fc-^p€p(s)g-ifc-Ao^e(*) _ Q^k -r ^ik X^i^is) ^-ik-\pi^{t) \ ^ (147) 



This expression is similar to pU5() and can be calculated at low temperature in the same 
way (see appendix B). According to (|227|) . its leading term when k <C and /3 — > oo is 
also related to the polarisability Oh of the ground state of the hydrogen atom: 

-4^Pi'«H(^2^), (148) 

where we used ^^{k) ~ . Using (|in2|) . the third graph in (|143j) is easily seen to give a 
contribution identical to that of the cluster C(l, 1), with an additional factor — (k) / {k'^ + 
K^{k)) due to the "dressing" by the root cluster Cf, = C(1,0) or C(0, 1). The last graph 
in p4l'{j) involves similarly two such dressing factors. We find therefore, at leading order 
and for /cA <C 1, 



Xii(fc) ^ -47rpif«^ 1 _ -f—. + ( __ ) _4Vd^^( ___ ) _ (149) 



2^2 / -K^ x2 



k"^ + V ^2 -)- ^2 



k"^ + K? 



This expression vanishes at A: = thanks to the screening factors k"^ /{k"^ + k^). Adding 
H149|) to ()1U3() . we obtain the following approximation to x(^) ™- the atomic limit 

^2 ,^2.2 

x(fe)-- ,2 , 2 -47rp>H ,0 , 2 ' A;A««1. (150) 



Eq. ()15U() describes simultaneously the screening due to free electrons and protons (first 
term), and the screening due to the polarisability of the hydrogen atoms. A plot of (|15()|) 
is shown on Fig. [3 

The approximation (|150|) can be compared with the extended RPA-dielectric func- 
tion introduced by Ropke and Der ^Sj. We first note that, in the Maxwell-Boltzmann 
approximation, the RPA-dielectric function 

a=e,p ^ ' -J rUp ^p-k 

-Ep*' = fi^p^ /{2ma), is simply related to the function K?{k). If we replace the factor (|4U() 
in (|56j) by (|39j) and perform the "time" integral, we find indeed eRPA(fc) = K^{k)/k'^. 
Using ®, eRPA(^) is therefore equivalent to the mean- field response function H1U3|) : 

^ l = XMY{k). (152) 



eRPA(fc) 



^The part (j)^'^'^'^\Ca, Ch) associated to the non zero frequency terms in (|93(l is integrable at zero density 
(see 20]), so that its contribution to p45l) is 0{z) for all values of k (discarding powers of (3). 
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Figure 2: Plot of the response function x{^)j according to H150|l . The electron-proton plasma is in an 
dilute atomic phase of density pif — 10^^ atoms/cm'^. The temperature is T = 1350 K, and the Saha degree of 
ionization 7 = pjf /(pjf + Pa'^) = 10^^^. In this situation, the cross-over distance l|47(l between ionic and dielectric 
screening takes the value Ai ~ 700 cm. On the plot, the inverse Debye screening length k ~ 10^^ cm^^ is almost 
at the origin, and the inverse lengths l/(/3e^) ^ A^^^ ~ 10® cm"-'^ are far to the right. This plot of xC^) shows 
the plateau of this function for k in the range X^^ <C A; <C A^^^, which corresponds to the dielectric response of 
a gas of atomic hydrogen. 



The extended RPA-dielectric function includes a contribution eat(fc) associated to the 
polarisability of the hydrogen atom: 

e(fc) = l + eRPA(fc)+eat(fc) (153) 

with 

n,n'=0 
n-n'=0 

{En, Ep and |j4nn'(^)P defined in section IH}. The above expression of eat(^) follows 
from eq. (2.4) of ^H] (see also eq. (4.265) in ^S]) by setting lu = 0, neglecting degeneracy 
effects, and retaining only the terms involving at least one ground state (because e"'^^" ^> 
e~^^", n 7^ 0). Comparing (|154() with H124() and H39|) . we see that eat(fc) is the same as our 
expression xit^(^) l|125|) . The approximations (|15flj) and (|153j) for the response function 
are therefore equivalent at leading order: they differ only by a term oc (p^an)^- Ropke 
and Der derived eat(fc) in the framework of the Feynman perturbation theory of the many- 
body problem. They must hence sum an infinite number of graphs to obtain the atomic 
contribution, since chemical binding is a non-linear effect in the interaction. This is in 
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marked contrast with our non-perturbative approach, where the atomic contribution (and 
the contributions of other bound states) are described by simple graphs. 

We stress that the approximation (|150j) does not include in a systematic way the 
leading corrections to XMF{k) for all values of A: <C Xi^, even in an atomic phase satisfying 
the inequalities (|21j) . Consider for example the contribution of the graph made of a root 
cluster C(l,2) describing a ion. According to H129I) and ()130|) . this graph will lead, 
when suitably dressed, to a term of order 



When k K, (|155|) is hence of the order of a polynomial in /? times Pj^/P^^ (since cc z). 
Although Pyz/ P^^ <C 1 in an atomic phase satisfying (|21() . it is much larger than p'^au 
when /? oo, as can be seen by inserting the values of the binding energies Eg, ~ —13.6 eV 
and £^12 — —14.3 eV. The second term in the approximation H15U|) is therefore not the 
leading correction to x(A;) — xmf(^) when A: ~ k. Obtaining a coherent approximation for 
x{k) at low densities for all values of k is thus a subtle problem, since all ions contribute 
to the constitution of the screening cloud when k is in the cross-over region. 



In conclusion, the screened cluster expansion is a convenient technique to study at low 
densities the response function (and other static equilibrium quantities) of the quantum 
plasma when the charges recombine into bound entities. With this technique, we have 
been able to determine the wave length range where x(fc) shows a plateau corresponding to 
dielectric behaviour. For k in this range, we have calculated the dielectric constant up to 
first order in the atomic density, taking into account in a controlled way all effects induced 
by the Coulomb potential (chemical binding, collective screening, polarization). The 
present method allows in principle quantitative calculations of other type of contributions, 
coming from higher order density terms, thermal excitations, more complex chemical 
species, etc. Finally, it can be extended to the study of more general partially ionized 
nucleo-electronic plasma. 



A Upper bounds on remainders 

In this appendix, we prove that the various terms we have discarded in section 6 are indeed 
negligible in the atomic limit. The points (i) and (ii) are established in section A.l and 
A. 2 respectively. Section A. 3 contains a few lemmas which are used in the proofs. The 
letter C denotes throughout this appendix a positive constant which can have different 
values at different places. 



A.l Neglecting screening effects in the atomic contribution 

We seek for upper bounds on the contributions to x{k) of the terms Ri (see ()lU9p - Hlll() ). 
which arise when the screened potential is replaced by the bare Coulomb potential in (|in5|l . 
These contributions have the generic form 



47r/3e2 



dr / D(C) / D(^ 



ds- 



2zp 



2z„ 



/o (27rA2)3/2 (27rA2)3/2 

^ gifc Ap4p(.)(gifc.[r+Ae4e(s)-Ap$p(s)] _ l)(e-ifc-^ _ l)R{r,Xe,Xp), (156) 
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where R = Ri, R2, Rs- We introduce absolute values, use |e*^ — 1| < C\x\ and | exp[— /3I/j(£)]| < 
C (this follows from 1)87(1 and (|97() ). Expression (|156() can be majorized by 

"^"'(2^(2^' I * / °(«»' / °«") /' ^'^^ 

X |r + A„C(«)-Ap4p(s)||r|. (157) 

We study first the contributions of the terms R2 and R3 , and then that of i?i . 



Contribution of R2 

Inserting I|11U() in ()157() , we have to estimate a finite sum of terms of the form 

P{P)z,Zpl dr I D{U I ds\4>,^{r,Xc,Xpr\V{r,Xe,Xp)r 

X |r + Ae^e(s)- Ap^p(s)||r| (158) 

where -P(/3) is a polynomial in /? and n + m > 6. It is proven in ^O" that, at low density, 
the potential (pch can be majorized by a constant times the electrostatic potential (|66() 

\(pchir, Xa,Xb)\ <C\Vciec{r,Xa,Xb)\- (159) 

Using this majoration, and introducing the change of variables r = \x, we obtain the 
following upper bound for H158|) : 

P{^P)z,z^j dx\x\ I D(0 / D(^p)(^'ds|a^+y^C(s)-y^^p(5) 

dTa / dn- 



^ + yS^e(ra)-y^^p(r, 



dr- 



1 



I^ + a/^Ut) 



(160) 



Let us show that the above integrals are finite, thereby proving that the contribution of R2 
is dominated by a polynomial in times ZgZp. 

Using twice the Cauchy-Schwartz inequality on the measure D(^g)D(^p), one can ma- 
jorize the integrals in ()160j) by J dxJ{x), where 



J{x)^\x\J ds\x + ^{s)\)\ 



4m 



X / D(C) / D(^ 



dTa / dn 



1 



An 



■ (161) 



^ + V^^e(ra)-^^^p(r,) 

In (jl61j) . we have used to write the functional integrals involving a single time in 

terms of the relative path ^(s) = m/me^g(s) — y^m/mp^p(s). The following lemma can 
be used to prove the finiteness of the integral over J{x). 

Lemma 1 (smoothing the Coulomb singularity) 



ds 



1 



ma) / D(e 



x + i{s) 



dT„ 



< 



C 



\xr + 1 



dn- 



< 



(162) 

c 

(163) 
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The functional integrations over the Brownian paths smooth out the Coulomb singularity 
at the origin so that H162() and (|163|) remain finite when \x\ 0. Note furthermore that 

J B{^)(^l\s\x + ^{s)\y <C\x\\ (164) 

The function J{x) can therefore be majorized by J{x) < C\x\'^ / {\x\'^^"'~^"^^ + l)"*^^^, and 
J dxJ{x) is hence integrable since n + m > 6. The contribution of R2 to x(^) is thus 
indeed dominated by a polynomial in ^/]3 times z^Zp. 

Contribution of 

We seek for an upper bound for 

PiP)z,Zpl dr I D(0 I D(^p)^'ds|r + Ae^e(^)-Ap^p(s)||r| 

00 ^ 5 
X i;-(/3|</'ch(r-,Xe,Xp)l)" Y.(P\V{r,Xe,Xp)\r. (165) 



n=6 m=0 



We write |(/'ch|" = |</'ch|" ^|</'ch|^ and, using (|97jl. majorize the sum over n by 
00 ^ 

E-(/3|'^ch(r,Xc,Xp)l)"-'<e^l'^^'^l <C. (166) 



After this majoration, we are left with an integral already studied (see HISHI) ). so that the 
contribution of R3, like that of R2, is dominated by a polynomial in ^/]3 times z^Zp. 

Contribution of Ri 

Since |-Ri(r, Xe, Xp)l can be majorized by C/3e^K;Sy^(r, Xe, Xp)) using H1U9|) . (|U7|) and the 
positivity of ;Sy'^(r, Xe, Xp) (see pi7|) ). we expect its contribution to x(^) to be smaller by 
a factor /3e^K as compared to the leading contribution Xat(^) (!116j) . Writing (|157j) in the 
center of mass and relative coordinates, we find that the contribution of Ri is dominated 
by C/3e2K(27rA2)-3/2exp[2/3;u] times 



(27rA2)3/2 



JdrJ D{^) l\s\r + Xm\\r\Bl.l{r,X^). (167) 



We show here that (|167j) behaves like a polynomial in times exp [-PEo], just as ^lE^. 
The contribution of Ri is hence dominated by P{\/p) {Pe"^ k) p^^, which is exponentially 
smaller than Xat(fc)- 

In order to prove that H167() behaves like exp[— /3i?o] at low temperatures, one must 
convert the functional integral into operator's language to extract the ground state en- 
ergy Eq . The main difficulty in the proof is then to keep track of the convergence of the 
r-integral, which is obvious in the functional integral representation, but is non trivial in 
the operatorial expressions. In the language of operators, p67|) is equivalent to 



e2 



drTr{u(/3-r)|q|U(r)|q|-e-^^^oT[|q|(r)|q| / dsV(s)) ]}, (168) 

where we use the same notation as in (|119|) . From the Dyson series, we have 

\J(t2 - ti) = e-(*2-*i)(H+V) = exp[-i2Ho]Texp [ - / ' ds V(s)j exp[tiHo], (169) 
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so that (|l(i8j) can be rewritten as 



IqlMlqlfe-Zo'^'^^VW.^i. 



n=0 



dsY{s] 



(170) 



The truncation terms in B'y2 therefore subtract out of the trace the terms of order 
to 5 in V, ensuring thereby the finiteness of the trace (the same argument shows that 
the trace pi9|) converges). To exploit these cancellations, we introduce in H168|l a limited 
Dyson expansion of U(s) to fifth order in V: 



U(S) = ^UW(5) + U(°)(5), 

fc=0 

where U*^'')(s) = Uo(s) = e"^^'' is the free propagator, 

U«(s) = Uo(s)-^t( ds' V(s'))', A; = 0, 1, . . . , 5, 
with V = —V = e^/|q|, and the remainder is 



(171) 



U(D)(^ 



dS2 



dS6 Uo(s)V(si) . . . V(S5)U0(S6)VU(S6). 



(172) 



(173) 



Notice that the \J^^\s) involve only free propagators Uo(s) (and Coulomb operators V), 
while U^^^(s) contains the full evolution operator U(s). Because of the truncation, the 
terms in ()168() of order less than 6 in F cancel out, and we are left with terms of the form 



^%rTv{u«(/3-r)|q|U('")(r)|q|}, 



(174) 



where Z, m = D or 0, 1, 2, 3, 4, 5 with 6 < / + m < 10 if /, m 7^ D. 

Consider first the case where / and m are different from D. Since only free propagators 
Uq are present, no bound state can occur in these terms. Let us show that ()174|) is then 
bounded by a polynomial in implying a bound P{^/0)zeZp for the contribution of 
these terms to x(fc). Inserting ()172l) into H174|) . we find 

-— I dr dsi... dsi / ds[ . . . ds'^ 

Tv{Uo(/3)t[v(si) . . . V(sO|q|(r)V(4) . . . V(4)|q|] }, (175) 

which can be majorized by 

^ dr dsi... dsi+m Tr{Uo(/3)T [|q| (r)|q| V(si) . . . V(sz+^)] } , (176) 

by extending all the integrals to the domain [0, /3] (the integrand is positive, as can be seen 
from the functional integral representation of the trace). Expression (|176() is equivalent 
to 



2\l+m+l 



1 



l\m\ (2^A2)3/2 



dr / D(^) / dr|r + A^(r)| \r\ 



ds 



1 



r + A^(s) 



l+m 



(177) 



By the same argument (based on lemma ^ already used in estimating the contribution 
of i?2) we conclude that H177|l is bounded by a polynomial in y^. 
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Consider now the terms with at least one index equal to D. We use the invariance of 
the trace under cychc permutations and the fact that |TrA| < ||A||i (where || • ||i is the 
trace norm), to majorize (|174|) by 



e2 



%r|||q|U«(/5-r)|q|U(-)(r)||i. (178) 







One can obtain the fohowing results on the norms of the operators 
Lemma 2 ( Convergence of the traces ) 

a) |||q|U«(s)|| <P(V^), / = 1,2,...,5 (179) 
h) |||q|U(°)(s)||i <P(^/^)e-^^o (180) 
c) |||q|Uo(t)|q|U(°)(s)||i<P(Vi)P(Vt)e-'^^». (181) 

The operator norm in lemmata holds only when I > 1, because U('-°)(s) does not 
contain a Coulomb operator to counterbalance the presence of the unbounded opera- 
tor |q|. Naively speaking, the trace norms in|2jb and|2lc are finite because U(°)(s) is of 
order 6 in V, which is enough to ensure convergence in 3-dimensional space even with two 
operators |q|. Furthermore, the bounds [Sb and[2jc grow like exp [sEo] because U(°)(s) 
contains the full evolution operator U(s) which involves bounds states. 

From the upper bounds of lemma and the inequality ||ST|| < ||S|| • ||T||i, one easily 
deduces that (|178|) is dominated by a polynomial in y/j3 times exp[— /?£'o], as announced 
after H167|) . Indeed, if one index is equal to D (say m = D), one uses ()181() if Z = 0, 
and (|179|) together with (|18()|) if / = 1,2,..., 5. We use also of course the fact that 
P — T and T do not exceed /? and majorize P(\//3 — t) and P{^/t) by polynomials P^^fjS). 
The same argument applies to the case / = D, m 7^ D by using the invariance under 
cyclic permutations in the trace (|174() to bring U*^^^ (/? — r) to the right of the expression. 
Eventually, the case / = m = D is treated similarly by applying twice (|18fl|l together with 
||T|| < ||T||i. 

A. 2 Neglecting excited and ionized states contributions 

From its definition 1)1221) . the contribution of excited and ionized states to x(fc) is 
xil^(fe) = -Ane^e^f''' dr /,/^(fc)Tr{u(/3 - T)QB_fcU(T)QBfc 

B_fc(r)Bfc5]-(y^ dsV(s))"]} (182) 



- e-^«oT 



with 



Bfc = (e^'^ 'irq - e-''=-3r^)^ = e-''^'^'^ (e^'^''^ - 1) A- (183) 

\k\ \k\ 

Expression 1)182(1 involves a truncated trace similar to l|168() . except for the additional 
operators Q that project the wavefunctions on the set of excited and ionized states. Notice 
that the operator B^ plays here a role similar to |q|, because we look for an upper bound 
that is uniform in k and lim^^o Bfc = ik ■ q^. The methods used to derive upper bounds 
on (|168|) can also be applied here with minor modifications (see below). The result is that 
xit^(^) is dominated by P{\fi3)p'-^e~^^^'^~^''\ where Ei is the energy of the first excited 
state of the hydrog en atom. ^|^^^(fc) is hence exponentially smaller than Xat^(^)- 

Following the same method as in section A.l, we introduce in (|182|) the limited Dyson 
expansion ()171|) of U(s). We distinguish two types of terms: 
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• Class I: terms with at least one 

• Class II: terms without \J^^\s) (th ese terms are of order to 10 in V). 

Let us show that the terms in the class I are dominated by P{^/P){2ttXI)-''^/'^ exp[-/3(^i - 
2/i)] and hence are o(p!f). Majorizing the traces by trace norms (and using the cyclicity 
of the trace), these terms are dominated by 

^(2^^Alj^ ■ I "^'^^^'^^^ " ^)QB-fcU(™)(r)Q||i, (184) 
where at least / or m is equal to D (the other index running from to 5). 
Lemma 3 The following bounds holds uniformly in k: 

a) ||BfcU«(s)Q|| < ||BfcU(')(s)|| <P(Vi), / = 1,2,...,5 (185) 

b) ||BfeU(D)(s)Q||i <P(V^)e-^^i (186) 

c) ||BfcU(o)(t)QB_fcU(D)(s)Q||i <P(V^)P(Vt)e-'^^i. (187) 

As opposed to the bounds of lemma HI the above bounds involve the first excited state 
energy Ei of the hydrogen atom because of the operator Q. Applying on (|184)) these 
bounds in the same way as we applied the bounds of lemma [21 on ()178|) , we obtain the 
announced result (|184|) = o{p'^). 

We show now that the terms in the class II grow at most like P{y/]3) exp[2/3/i] and hence 
are also o(/3^). Contrary to the case considered in section A.l, the terms in class II of 
order less than 6 in V do not cancel out, because of the presence of the two operators Q. 
We write therefore Q = H — P and look first at the terms without any P. Just as in 
the contribution of exact cancellations occur now at order to 5 in V. The terms 
of order 6 to 10 in V (and still without P) are most easily majorized in the functional 
integral representation. Since | exp[— ifc • ^r](exp[zfc • r\ — l)/|fc|| < C\r\ uniformly in /c, 
these terms are dominated, in absolute values, by C^exp[2/3/x] times the expression (|177() 
which grows polynomially with y/]3. It remains thus only to consider the terms in the 
class II with at least one ground state projector P. Consider for example the term with 
\]^^\f3 — r)P (the other terms can be treated similarly). Using the cyclicity of the trace, 
we can majorize this term by the matrix element 

^(^^^AJV^ /o d^l^0|PB-fcU(™)(r)BfcU«(/3-r)|0)|, (188) 

where l,m = 0,1,..., 5. Thanks to the fast (exponential) decay of the ground state 
wavefunction, the operator PB_fc (and its adjoint B^P) are bounded uniformly in k. If 
/ 7^ 0, the matrix element involves therefore a product of bounded operators (see lemmaOJa 
and note that ||U(™)(s)|| < as can be seen from the proof of lemma [21a), and grows 
thus at most like a polynomial in y/^. If / = 0, one commutes B^ with U^'^^ {[5 — r) and 
rewrite the matrix element as 

(0| PB_fcU("*)(r)u(0)(/3 - r)BfcP |0) + (0| PB_fcU('")(r)[Bfc, U(°)(/3 - r)] |0) . (189) 

bmce [Bfc,UW(/3 - r)] is bounded (lemma [Ha), one can apply again the same argument. 
The terms in the class II grow thus indeed at most like P{\/]5) exp[2/3;u]. 

A. 3 Proof of the lemmas 

Here, we set units such that = 1 and /m = 1, so V = — l/|q| and Uo(s) = exp[— sHq]. 
To prove lemmas [21 and [HI we need a few basic results on norms involving the free propa- 
gators and the Coulomb operator: 
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Lemma 4 The operator Uo(s)V is bounded for s > 0; 

||Uo(^)V|| <2(i= + 2) (190) 

Lemma 5 The following commutators are bounded: 

a) l|[Uo(s),q;.]|| <CV^ (191) 

b) ||[Uo(s),q^]V|| <C(1 + ^) (192) 

c) ||[[Uo(s),q^],q,]V|| <CV^(1 + ^/^) (193) 

Lemma 6 The operator Uo(s)VUo(t)V belongs to the Hilbert- Schmidt class for s,t > 0; 

||Uo(5)VUo(t)V||2 < ^=£= (194) 
^/s{s + t) 

Lemma 7 The bounds of lemma [3| remain valid when ( and q^) is replaced by the 
operator |q| = + + q| . 

Lemma 8 The following bounds hold uniformly in k: 

a) ||[Bfc,Uo(s)]|| < C^/^ (195) 

b) ||[Bfc,Uo(s)Q]|| <C(1 + ^) (196) 

c) ||[Bfc,Uo(s)]V|| <C(1 + ^/^) (197) 

d) II [Bfc, [B_fc, Uo(s)]]V|| < CV^il + V^) (198) 

The fact that the operator Uo(s)V is bounded (lemma 2)) can be traced back to the 
uncertainty principle of quantum mechanics. Indeed, if V is very large, the wavefunction is 
close to the origin and hence well localized. By the uncertainty principle, the momentum 
must be large, so that Uo(s) = exp[— sp^/2] is small. This makes the product Uo(s)V 
bounded if s > 0. Concerning lemmalHl recall that the Hilbert-Schmidt norm is defined as 
||A||2 = VTrA^A. Naively speaking, lemmaElholds because the trace TrA^A converges for 
the following two reasons: the product Uo(s)V remains finite at short distances (lemmallj 
and is integrable at large distances. In lemma Ela, the presence of a commutator is 
crucial to ensure the finiteness of the operator norm (when s = 0, this norm vanishes). In 
lemma inib, the commutator does not play such a crucial role, because q^V is a bounded 
operator (||q/xV|| < 1). As already mentioned, the bounds of lemma|Hlare essentially the 
same as the bounds of lemma El and d because limfc^oBfc ~ fc • q. Lemmas 01 El and IHl are 
established in We prove below the lemmas 2, 3 and 7, 8. 

Proof of lemma (21 

a) We introduce H172|) in 

II |q|U('=^(s)|| < r dsi . . . dsk II |q|Uo(s - si)V|| 

Jo Jo 

X ||Uo(si - S2)Y\\ • . . . • ||Uo(sfc_i - Sk)Y\\ ■ ||Uo(sfc)||. (199) 
Notice that the operator |q|Uo(s)V is bounded (this follows from lemma [71b): 

|||q|Uo(s)V|| <C(l + ^/i). (200) 
Using lemmal^and ||Uo(sfc)|| < 1, (|199)) involves integrals of the form 

/ ' ' dsfc ^ ^ = 2^ < 2^^, (201) 

Jo \/Sk-l — Sk 
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Hence only integrable singularities are present, and || |q|U'^^''(s)|| < P2k{Vs), where P2k{Vs) 
is a polynomial of order 2k in ^/s. 

b) From ()173() . one has 

II |q|U(°)(s)||i < r dsi . . . r ds6 II |q|Uo(s - si)V|| 
Jo Jo 

X||U0(S1-S2)VU0(S2-S3)VU0(S3-S4)VU0(S4-S5)V||1-||U0(S5-S6)V|H|U(S6)||. 

(202) 

The bound ()18U() follows then from the use of ||AB||i < ||A||2 • ||B||2 (the analogue of the 
Cauchy-Schwartz inequality for operators), lemmalUEl (|2()()|) and ||U(s6)|| < exp[— sE'o]. 
Note that the integrals are convergent. In particular 

fS rsi rS2 rS3 

/ dsi ds2 dss dg4 ^^ ^= = 2s\ (203) 

Jo Jo Jo Jo ^/{Sl - S2){S - S2){S3 - Si)[S2 - S4) 



c) The trace norm in H181|) is majorized by 

r dsi . . . r dsQ II |q|Uo(t)|q|Uo(s - Si)VUo(si - S2)V|| 
Jo Jo 

X ||Uo(s2 - S3)VUo(s3 - S4)VUo(s4 " S5)VUo(s5 " S6)V||l • ||U(s6)||. (204) 
Proceeding as in the proof of the point b), it remains only to show that 

r dsi r ds2 II |q|Uo(t)|q|Uo(s - si)VUo(si - S2)V|| < P{y/^)P{Vt). (205) 
Jo Jo 

To prove (|205j) . we commute repeatedly the operators |q| to the right, making |q|V = U 
appear, or the bounded operator |q|Uo(s)V (see (|2()()j) ). We have thus, without specifying 
the time arguments, 

|q|Uo|q|UoVUoV= [|q|,Uo]|q|UoVUoV + Uo|q|Uo|q|VUoV 

+ Uo[|q|,Uo]|q|VUoV + Uo[|q|, [|q| , Uo]]VUoV. (206) 

The bound (|2U5() follows then by applying the triangle inequality on (|206() and using lem- 
mas m and [71 



Proof of lemma |31 

The lemma El can be proven exactly in the same way as lemma 121 if we use the bounds of 
lemmalSlin place of the bounds of lemmajTl The factor exp[— si?i] comes from ||U(s6)Q|| < 
exp[— se-Ei] < exp[— si?i]. 

Proof of lemma [71 

a) We use the upper bound ||A|| < max(S'3;, S'y), where 5a; = supg, / dy\ {x\ A |y) |, Sy = 
supy / da;| {x\ A \y) \ and A = [Uo(t), |q|]: 

f 1 i^-i/P 

Sx = Sy = sup / do; — — — e 2r ||a;| - |y|| . (207) 
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We introduce the change of variables x = u + y, majorize Sy by inserting supy under the 
integral sign and use supy + — \y\\ = \u\. We obtain thus 

Sy<-, \:tt^ I due-^\u\ <C^/T, (208) 



(27rT)3/2 

from which lemma [Tj a follows. 

b) Since |q|V = H, one has 

||[Uo(s), |q|]V|| = ||[Uo(s),q'V]V|| < ||Uo(5)q'v2|| + ||qVUo(s)V|| (209) 

where we used the triangle inequality. The first norm is smaller or equal to one, and we 
majorize the second norm by 

3 3 

||q'VUo(s)V|| < ^ ||g2vUo(s)V|| < ^ llg^Vll • |k;.Uo(s)V|| < C(l + (210) 

The last inequality follows from ||g^V|| < 1 and lemma El b. 

c) We use [[A,B],B] = [B^, A] + 2[A,B]B and the triangle inequality to write 

||[[Uo(.),|q|],|q|]V|| <||[q^Uo(.)]V||+2||[Uo(.),|q|]|q|V||. (211) 
From lemma |7|a, the second norm in ()211() is dominated by C^/s. Since 

[q': Uo]V = ^(2[q^, UoKV + [q^, [q^, Uq]]), (212) 

we deduce from lemma El a, El c and ||q^V|| < 1 that the first norm in (|211j) is bounded by 
C^/s{l + y^), and hence lemma |7|c is proven. 

Proof of lemma jH 

We give the main steps in proving the points a) and b). The proof of the points c) and 

d) are left as an exercise to the reader. 

a) It is enough to show that 

gifc q _ 2 



„. ,Uo(t)]|| <C^, Vfc. (213) 

The operator exp[zfc • q] performs a translation in the space of momentum. Introducing 
the new notation \J^\p) = exp[— Tp^/2] for the free evolution operator, we have 

UW (p)e*'=-'i = e*'^-'iU(.°) (p + k). (214) 

The commutator in (tTnT|l thus evaluates to (U^-°^(p-fc)-U^-°^(p))e*'' 'i/|fc|. The bound JMI) 
follows then from || exp[ik • q]|| = 1 and 

||u(°)(p-fc)-uW(p)|| ^ suppa3|e-^(P-'')'-e-iP^| ^ 

1^1 ~ 1^1 ~ 

In (|215j) . the last inequality is uniform with respect to k. 

b) Writing Q = II — P, we have from the triangle inequality 

||[Bfc,Uo(s)Q]|| < ||[i?fc,Uo(s)]|| + ||[i?fc,Uo(s)P]||. (216) 

The first norm is bounded by C^/s according to lemma |H|a. The second norm is ma- 
jorized by ||Uo(s)[fifc,P]|| + ||[Bfc,Uo(s)]P||. We obtain the bound (tTT^HIl by using ||P|| < 1, 
||Uo(s)|| < 1, lemmalHIa and ||[Bfc,P]|| < C. 
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B The dressed atomic contribution 



We calculate in this appendix the leading term in the atomic limit /? ^ co of the dressed 
atomic contribution (|147j) . by using the same method as in section l6?2l Replacing (j)hy V 
(using pU8() ). and omitting the factor —K?/{k'^ + k^), we have to evaluate 



Bvlir, Xe, Xp) (e^'=-^^^-(")e-^'=-^^^-(*) + e*'»-^p«p(^)e-*'=-^p^pW 

_ e~*'° '^e*'° '^p^p^''-'e~**' '^"^<=^*-' — e**'''^e**'''^'=^<=^'*-'e"*'° '^p^p^*-'y (217) 



We use the center of mass and relatives variables (|113j) . The integrations over the center 
of mass factors, and is easily calculated using the covariance (jSH: 

(2vrA2)3/2 J ^^^-^^ (27rA2)3/2^ - ^1--*^'^^ 

We hence obtain the expression 



4e*^ d«/„_„(.)^/ci./D(e) 



Xd.at.(fc) - ^2 

^vli'^, A^){e'^'=-^^«(^)e-^'=-^^^W + e-^'^-lf ^«We^'=-Tf 

_ Q-ik-r^-ik-^Xi{s)^-ik-!^Xm _ Qik-r^ik-^Xiis)^ik-^Xmy (219) 

The integrand is a symmetrical function of s and t, since ;Sy^(r, A^) = i3y^(— r, — A^) and 
D(^) = D(— ^). The time integrations can hence be limited to the region s > t: 

/ ds dt^2 ds dt. (220) 
Jo Jo Jo Jo 

It is convenient to factor the terms in the braces in 1)219^ . according to 

(Q^k-'-^lr+^s)] _ g-ifc.^[r+A4(s)] j j^^-ifc- ^ [r+A^{t)] _ ^ifc- ^ [r+A4(t)] j _ (221) 

In operatorial language, Xd.at. (^) becomes 

Xd.at(fc) = -^4e2/3^^ ds dt/|,_i|/^(fc)Tv{e-^HAL(^)Ai„t(t) 

- e-^"«T [At (s) A(t) ^—^[f d'uV('u)) "] } . (222) 

n=0 

where Aint(s) = U(,s)AU(— s) is the time evolved operator A (compare with (|119j) ). At low 
temperatures, the leading terms arise from terms with ground state contributions from 
the evolution operator U(r) = e""^^ (r = t, s,/3). Introducing the ground state projector 
P = |0) (0| and Q = H — P, the first term in the trace, i.e. e~'^^A[j^^(s)Aint(t), can be 
written as (using the cyclicity of the trace) 

U(/3 - s + t)(P + Q)AU(s - t)(P + Q)A. (223) 
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Retaining only terms with at least one ground state projector (as in (|12;-i|l ). we find 

xSt. (fc) ^ -^4e2/3/^| f\s-t\/(3ik) 1v{u(/3 - s + i)PAtu(5 - t)FA 

+ U(/3 - s + t)PAtu(s - t)QA + U(/3 - s + t)QAtu(s - t)PA}. (224) 

In terms of atomic eigenvalues and eigenstates, Xd^^t (^) becomes 

XZW = -4.p?|J| /''d<,£dte-^'='^2Ta-T)(|A„„(.)|^ 

+ ^ e-(-'«''"-'=«'|A„„(fc)|' + ^ e-(''-+"<'=--''»'|A„„,(fc)p) (225) 

m>l m>l 

where we factored out the atomic density (|13|) and used |Ann'(fe)P = |An'n(fc)P- Let us 
evaluate (|225() for k <C A~^. The first exponential can be approximated by one at lowest 
order and we find, after having performed the time integrations, 

n=0 

We recover hence the polarisability aii{k,f3) ()43() . so that eventually 

xSt.W =^ -4vrpLV, A:Aa«l. (227) 
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